Torsion Balance Investigation of 

Casimir Effect 



A thesis 
submitted for the Degree of 

Doctor of Philosophy 

in 

The Faculty of Science 

Bangalore University 



by 

G RAJALAKSHMI 



Indian Institute of Astrophysics 
Bangalore 560 034, India 

May 2004 




Torsion Balance Investigation of the 

Casimir Effect 



A thesis 



OO 

o 
o 

>>■ 
00 



cr: 



> 

m 
oo 



in 
o 

00 

o 



X 



submitted for the Degree of 

Doctor of Philosophy 



in 



The Faculty of Science 

Bangalore University 



by 

G RAJALAKSHMI 




Indian Institute of Astrophysics 
Bangalore 560 034, India 

May 2004 



Declaration 



I hereby declare that the matter contained in this thesis is the result of the investi- 
gations carried out by me at the Indian Institute of Astrophysics, Bangalore, under 
the supervision of Prof. R. Cowsik. This work has not been submitted for the award 
of any degree, diploma, associateship, fellowship etc. of any university or institute. 



Prof. R. Cowsik G Rajalakshmi 

(Thesis Supervisor) (Ph.D. Candidate) 

Indian Institute of Astrophysics 
Bangalore 560 034, India 
May 2004 



Certificate 



This is to certify that the thesis entitled ' Torsion Balance Investigation of the Casimir 
Effect' submitted to the Bangalore University by Ms. G. Rajalakshmi for the award 
of the degree of Doctor of Philosophy in the faculty of Science, is based on the results 
of the investigations carried out by her under my supervision and guidance, at the 
Indian Institute of Astrophysics, Bangalore. This thesis has not been submitted for 
the award of any degree, diploma, associateship, fellowship etc. of any university or 
institute. 



Indian Institute of Astrophysics, 

Bangalore 

May, 2004 



Prof. R Cowsik 
(Thesis Supervisor) 



Acknowledgements 



I thank Prof. R. Cowsik for his guidance and continuous support as my thesis advisor 
during the course of this work. I have learnt many aspects of research from discussion 
with him. 

I thank Dr. C. S. Unnikrishnan for his support and encouragement throughout - as 
a colleague and as a friend. His never ending stream of crazy-sounding ideas are a 
constant source of inspiration. 

I have gained a lot from discussions with Dr. N. Krishnan and Prof. S. N Tandon. 
Their open criticisms have put me on the right track many a times. 

I thank Prof. R. Srinivasan, for the technical help and support he provided in building 
the experiment, especially the CCD optical lever. 

I thank Dr. B. R. Prasad and Dr. Pijush Bhattacharjee, who as doctoral committee 
members helped in advancing the thesis. 

I thank Prof. B. P. Das, Dr. Sharath Ananthamurthy, Prof. C. Sivaram, Dr. A. K. 
Pati and Dr. Andal Narayanan for the interesting discussions I have had with them. 

I thank Surcsh for his patient support and companionship during the years it took to 
build the experiments. 

I thank Ms. A. Vagiswari, Ms. Christina Birdie and all the other library staff for the 
excellent library facility provided at the Institute. 

I thank the present and past members of the Board of Graduate Studies for their help 
and guidance. 

I thank Mr. J. P. A. Samson, Mr. Thimmiah, Mr. Sagayanathan and Mr. Periyanayagam 
of the mechanical workshop for their help with the design and fabrication of the var- 
ious mechanical parts of the experimental apparatus. 

Building a lab requires support and advice from many people on aspects ranging 
from instrumentation to administration, it is impossible to acknowledge everybody's 



contribution individually. I thank all the staff of IIA who helped in the setting up of 
the 'Zerolab'. I thank the staff of the computer center, electrical section, photonics 
and the administration for their continuous support. 

I thank the present and past Chairmen of the Physics Department of Bangalore 
University and the members of staff in the University for their help with the admin- 
istrative matters related to the thesis. 

The several years I have spend in IIA was made pleasant by the many wonderful 
friends that I was lucky to have. I thank - Ramesh, Krishna, Sankar, Sivarani, 
Geetha, Sridharan and Rajesh Kumapuaram for introducing me into IIA; Mogna 
and Holger for being there during hard times; Arun, Rana and Pratho for those 
numerous discussion - scientific and otherwise; Dharam for being a nag and for the 
care; Manoj for the enumerable arguments on wide ranging topics; Appu for being 
such a patient listener; Ravi for help in a many a things- personal and technical; 
Preeti, Latha, Sivarani and Shalima for those many early morning chais; Geetanjali 
for being a wonderful roommate. I also thank E Reddy, Sujan, Annapurani, Srikanth, 
Rajguru, Pandey, Swara, Bhargavi, Sanjoy, Pavan, Mangala, Jana, Ramachandra, 
Ambika, Shanmugham, Kathiravan, Reji, Sahu, Malai, Jayendra, Nagaraj and Maiti 
for making my stay at IIA enjoyable. 

I thank Pramila, Rani and Samson for their hospitality and care. 

I thank my loving parents whose guidance and encouragement has made me what 
I am. My brother, Sankar for his silent support and affection. I thank Raja for 
enduring me through these years and being a constant source of encouragement and 
companionship. Athai and Athember for their affection and support. Last but not 
the least, I thank my grand fathers - Chandrasekharan, Ananthanarayanan and Kr- 
ishnamurthy - who at a very young age guided me in developing a personality of my 
own. To them I dedicate this thesis. 



CONTENTS 



1. Introduction 

1.1 Casimir Force - an introduction 

1.1.1 Casimir force as a manifestation of zero point energy 

1.1.2 Effect of Finite Temperature 

1.1.3 Effect of Finite Conductivity 

1.1.4 Importance of these effects for experiments 

1.2 Historical Review of Experimental Status 

1.3 Motivations to study Casimir force 

1.3.1 Understanding the Quantum Vacuum 

1.3.2 "The Hierarchy Problem" 

1.3.3 Constraints on new macroscopic forces 

2. Torsion balance- Design and Fabrication 

2.1 General principle of the apparatus 

2.2 Torsion Pendulum 

2.2.1 Mass Element 

2.2.2 The Pendulum Suspension 

2.3 The Capacitor Plates and Torsion Mode Damping 



2.4 The Spherical Lens and the Compensating Plate 

2.5 The Vacuum Chamber 

2.6 The Electrical Wiring and Grounding in the Apparatus 

2.7 The Thermal Panels 

3. The Autocollimating Optical Lever 

3.1 Conceptual aspects of the design 

3.2 Construction of the Optical Lever 

3.3 Implementation of the Centroiding Algorithm 

3.4 Tests and characteristics of the optical lever 

4. Characterization of the Apparatus 

4.1 The Torsional Pendulum 

4.2 Pendulum behaviour during evacuation 

4.3 Calibration of the force arm 

5. Casimir Force Measurement - Strategy, Data Acquisition and Analysis 

5.1 Strategy of the experiment 

5.1.1 Summary of Forces acting on the Pendulum 

5.1.2 The Effect of these Forces on the Pendulum 

5.2 Data Acquisition: 

5.3 Analysis: 

5.4 Discussion of results 

6. Comparison of experiments on Casimir force 

6.1 Experiments to Study Casimir force 



Contents 



V 



6.2 Comparison of results 

7. Bounds on the Strength of New Macroscopic Forces and Future Direc- 



92 



Appendix 

A. Casimir Force between infinite Parallel Plates 

A.l Force between Dielectrics 

A. 2 Casimir force at finite temperature .... 



B. Casimir Force between spherical lens and flat plate 

C. Yukawa Force on a Gold Coated Lens due to Gold Coated Plate . 



96 



96 



tions 

7.1 Constraints on new macroscopic forces 

7.1.1 Astrophysical Bounds 197 

7.1.2 Bounds from Laboratory tests of Inverse square law |9 

7.1.3 Bounds from Casimir Force Measurements 

7.1.4 Estimation of bound from our experiment 

7.2 Future Directions 1105 



99 



100 



109 



110 



114 



121 



124 



127 



1. INTRODUCTION 



Abstract: This chapter presents an overview of the theoretical background and the mo- 
tivations for the experiment. The chapter begins with a general description of Casimir 
force, with discussions on the effect of finite temperature and finite conductivity. This 
is followed by a short historical review of the earlier experiments to measure Casimir 
force. The recent motivations to study Casimir force are then presented. 

1.1 Casimir Force - an introduction 

Even Schwinger, who was awarded the Nobel Prize for his pioneering contribution to 
Quantum Electrodynamics, has remarked that one of the least intuitive consequences 
of Quantum electrodynamics is the existence of a force of attraction between two 
perfectly conducting uncharged plates |l5]. In 1948, H. B. Casimir [3], showed that 
for two infinite parallel plates, separated by a distance d, this force per unit area is 
given by. 



The force is independent of the charge or mass of the plate. For plates of 1 cm^ area 
separated by 1 fim, the force is comparable to the gravitational attraction of two 
400 g masses separated by 1 cm or the Coloumb force on the electron in a hydrogen 
atom. 



Fc{d) 



240^4 
0.013 



dyn. cm ^ where c?^ = c? in microns 



(1.1) 
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£^ z = d 

Fig. 1.1: Two perfectly conducting infinite plates, placed distance d apart. 
1.1.1 Casimir force as a manifestation of zero point energy 

When electromagnetic field is described quantum mechanically, it has properties sim- 
ilar to an assembly of quantized harmonic oscillators. Each mode of the electromag- 
netic field defined by a set of parameters like frequency and polarization, is represented 
by one oscillator. The allowed energy levels of an electromagnetic wave of angular 
frequency C(j are given by the Planck relation = {n + ^)huj,n = 0,1,2,3.... The 
integer n, for the electromagnetic field, corresponds to the number of photons. Thus, 
even the zero quanta or the 'vacuum' state of the field still contains field fluctuations 
which result in a non-zero energy. In free space, all modes (frequencies) of the elec- 
tromagnetic wave are possible and each mode has a finite energy. As a consequence 
the vacuum or the zero photon state of the electromagnetic field in free space has 
infinite energy and infinite energy density. However, physically real effects arising 
from quantum fluctuations in the vacuum turn out to be finite and "renormalized" . 
One example is the Lamb shift of atomic energy levels. 

These vacuum field fluctuations also give rise to measurable mechanical effects on 
macroscopic systems which manifest as Casimir forces [25] . If an infinite, planar cav- 
ity comprising of perfect conducting plates separated by a distance, d, is placed in 
'vacuum', it imposes boundary conditions on the zero-point electromagnetic fluctua- 
tions. As a result, the possible modes of the electromagnetic field is restricted within 
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the cavity. Thus there is a finite difference in the energy of the vacuum field outside 
and inside the cavity. This results in a quantum vacuum pressure that attracts the 
cavity plates together. For the zero photon state, this pressure is given by Eqn. 11.11 
(see Appendix A for details). 

1.1.2 Effect of Finite Temperature 

The non-vacuum state of the electromagnetic field also has a Casimir force associated 
with it. In general at any finite temperature, T, thermal fields are also present in 
addition to the vacuum electromagnetic fields and the possible energy levels are given 
by the Planck's spectrum. 



From these results note that the distance dependence changes from 1/d'^ at low tem- 
peratures to 1/d^ at high temperatures. 

Considering Eqn. 11.51 it is obvious that the important non-dimensional parameter, 

that distinguishes the domains of high and low temperature, is x = ATckBTd/hc. 

Moreover, it is interesting to note that high temperature also corresponds to larger 

^ The prime over the summation symbol means that a factor half should be inserted for the n = 
term. See for example [S] 





The Casimir vacuum pressure (force per unit area) is then given by [5J 




with C(3) = 1.20206 
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Fig. 1.2: Two semi-infinite dielectric slabs of dielectric constants ei and 62, placed 



separation d between the plates and vice versa. Thus at any given temperature, the 
law which governs the vacuum pressure will depend on the distance between the 
plates. No experiment till date has been able to observe these finite temperature 

corrections to Casimir force, due to limitations in sensitivity in the distance range 
where such effects start to become significant. 

The primary aim of the work described in this thesis is to observe the Casimir force in 
the distance range where this change in the distance depcndance of the force occurs. 
At a temperature of ~ 300° K, the change from l/d^ to 1/d^ is expected to occur at 
about 2 //m to 4 //m. Our experiments scan separations from about 1 //m to 10 //m 
and thus will be able to probe this change over from the low temperature to the high 
temperature domain. 

1.1.3 Effect of Finite Conductivity 

The discussions so far assumed that the cavity plates are perfectly conducting, i.e, 
they have infinite conductivity at all frequencies of the electromagnetic field. In 
experimental situations, this simplified assumption is unrealistic and the dielectric 
properties of the cavity plates should also be considered. 



distance d apart in vacuum. 
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Liftshitz [1] , developed the first macroscopic theory of forces between dielectrics. His 
results reduce to the Casimir force given by Eqn. Il.ll for the case of perfect conductors. 
For finite size plates with dielectric constants €1(00) and €2(00) [Fig. II. 2j . the Casimir 
force per unit area is given by the Liftshitz formula, 



{si-p S2-p ) 

where Si, p are variables that depend on the dielectric constants of the medium 
and the wave vector k of the electromagnetic field, (see Appendix A for details). 

1.1.4 Importance of these effects for experiments 

The experiments on Casimir force are typically carried out at room temperature us- 
ing metals with finite electrical conductivity. In order to make a comparison between 
the experiments and the theory, it is essential to quantify the effects of temperature 
and conductivity. As summarized above the effects of the thermal field fluctuations 
on the Casimir force are known to become important when the spacing between the 
boundaries is of the order of the characteristic length, = cut is the dominant 
thermal angular frequency at the temperature T. Similarly, the plasma frequency Up 
of the metal determines the length scale, Xp at which the finite conductivity effects 
are appreciable, i.e, Xp = —. 

These effects have been calculated by several methods (see for example, [El [T71 17]). 
Fig. 15. 21 represents one such calculation [l7j. The figure shows a plot of the correction 




(1.7) 
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Fig. 1.3: Force correction factors as defined by Eqns. 11.81 11.91 and II. lOl are plotted as 
a function of spacing between plates at a temperature of 300° K (from, [IT]). 

factors estimated using, 

T 

■q 

V 

Thus 7]^ and rj^ are the estimates of the thermal and the conductivity effects respec- 
tively, rj is the combined corrections. As is clear from Fig. 15. 2[ the thermal effects 
and finite reflectivity of the metals are significant at quite different distance scales. 
At the sub-micron distances, the finite temperature effects are negligible, while the 
conductivity effects are considerable. Above a few microns the reflectivity effects are 
small as compared to the appreciable thermal effects. Thus, to study the effects due 
to finite temperature on the Casimir force, experiments need to be performed in the 
distance range of a few microns. 

1.2 Historical Review of Experimental Status 

Experiments to detect forces between flat plates have been attempted since the 17*^ 
century [17]. These early experiments were aimed at studying the existence or non- 




FM 

T p 

= r] 7]^ 



(1.8) 

(1.9) 
(1.10) 
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existence of the physical "vacuum". They studied the adhesion forces between two 
plates in 'evacuated' containers. In the 20*'* century, with the emergence of theories 
concerning long-range, London - van der Waals interactions and Casimir force, re- 
newed experiments were carried out. (For a review of the experiments until 2001, 
see 

The earliest attempt to measure Casimir force was by Overbeek and Sparnaay [15j in 
1952. They tried to measure the force between two parallel polished flat glass plates 
with a surface area of 1 cm^, in the distance range of 0.6 fim to 1.5 fim. The mea- 
surements at 1.2 fim, 'pointed to the existence of a force which was of the expected 
order of magnitude' [IT]. 

Derjaguin and Abrikossova [HI [7] were the first to obtain results in the distance range 
0.1 nm - 1.0 /im that were in agreement with Lifshitz's theory. Sparnaay [16] re- 
peated his measurements with metal plates in 1957. He measured the force between 
chromium plates and chromium-steel plates. The measurements did not 'contradict' 
the expected force per unit area from Casimir's relation. 

The next major set of improved measurements with metallic surfaces were performed 
by van Blokland and Overbeek [18] nearly 20 years later, in 1978. They measured the 
forces between a lens and a fiat plate coated with chromium using a spring balance 
at distances between 0.13 fim and 0.67 fim. This measurement can be considered as 
the first unambiguous demonstration of the Casimir force between metallic surfaces. 

In the last decade, attempts to understand the nature of quantum fluctuations at 
macroscopic scales and predictions of new long range forces at the sub-millimeter 
scales by theories that unify the fundamental forces have rekindled interest in Casimir 
force measurements. The earliest of these was by Lamoreaux in 1997 [24J, who mea- 
sured the Casimir force between a lens and a plate using a torsion balance in the range 
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0.6 fim - 6 fim. In a later experiment, Mohideen et al. [13] measured the Casimir 
force for separations from 0.1 fim to 0.9 fim using an atomic force microscope. Ex- 
periments at Bell-Labs by Chan et al. [1] indicate that Casimir type forces play an 
important role in micro-electromechanical systems. Recently Bressi et al. [2j have 
carried out high precision experiments between parallel plates in the range 0.5 /im - 
3 /im and the related force coefficient was determined at the 15% precision level. The 
most recent experiment by Decaa et. al. measures the Casimir force between two 
dissimilar metals for separations of 0.2 /xm - 2 /^m. 

A more detailed review and a discussion of data from the recent experiments on 
Casimir force will be presented in a later chapter. All the experiments performed so 
far to measure Casimir force were carried out at room temperatures and probed the 
distance range of 0.1 fim to 3 /im. As the separation increases, the force decreases 
rapidly as d~'^ to start with and as in the finite temperature regime [Eqn. II. 5j . It 
is important to measure the force at (i > 3 fim to detect and characterize the finite 
temperature corrections to the Casimir force. 

1.3 Motivations to study Casimir force 

Historically Casimir derived his results while attempting to explain the inter-molecular 
interactions seen in experiments on colloidal suspensions. Casimir force was looked 
upon as the effect of finite speed of light (retardation effect) on the London - van 
der Waals interaction. The early experiments that measured the force of attraction 
between surfaces where aimed at understanding the inter atom interactions and to 
see the change over from the van der Waals force at very small separations to the 
retarded van der Waals or Casimir force as the separations increased. The interest 
slowly waned once the macroscopic theories of interactions were verified by experi- 
ments. 

In modern times, interest in Casimir force has been aroused by the crucial role it plays 
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in theories of fundamental physics. Casimir force provides exphcit evidence for the 
existence of vacuum fluctuations and for the interplay between the microscopic (quan- 
tum) and the macroscopic worlds. Also, in the last couple of decades, several new 
theories have been proposed that predict new physics in the sub-millimeter distance 
range. Casimir force is the dominant background in this range and hence it becomes 
essential to understand all aspects of Casimir force before looking for new forces at 
these scales. 

1.3.1 Understanding the Quantum Vacuum 

The existences of electromagnetic field fluctuations in vacuum presents problems due 
to the amount of energy it carries. General Relativity states that all forms of mat- 
ter and energy should gravitate. The large energy density, associated with the 
vacuum fluctuations should induce very large gravitational effects, much larger than 
that allowed by observations. This "vacuum catastrophe" is related to the famous 
cosmological constant problem (see [SB [S21 HH SO] for a review). 

In 1917, when Einstein first attempted to apply his new theory to relativity to the 
Universe, the Universe was believed to be static. Einstein could not construct a static 
universe if there was only matter and curvature, so he introduced a free parameter A, 
the cosmological constant, into his theory which was a form of energy with negative 
pressure. With the discovery of the expansion of the Universe, Einstein suggested 
that this could be dropped. But it was retained alive in discussions of cosmology, 
and has been used time and again to explain observations that did not fit into the 
standard scenarios in cosmology [IHl HE]- Currently, there is strong observational 
evidence W\\ for an accelerated expansion of the Universe. This can be explained 
by a non- vanishing cosmological constant. Thus, in the present scenario, the geometry 
of the Universe is determined by the energy density of matter, p^, the energy density 
due to vacuum, and that from the cosmological constant, A/(87rG'). The quantum 
vacuum has properties similar to those attributed to a positive cosmological constant. 
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the most important property being an effective repulsive gravity. This is because 
the acceleration is proportional, in General Relativity, to the term — (p + 3p) of 
the matter, where p is the energy density and p the pressure. For normal matter, 
the term (p + 3p) is positive and therefore, the universe is expected to decelerate 
as it expands. The diagonal elements of the energy momentum tensor are energy 
density and the three components of pressure. Therefore, an energy momentum 
tensor that is proportional to "vacuum" with diagonal elements (1,-1,— 1,-1) will 
have its equation of state p = —p, and the effective acceleration of the Universe 
with such a source term will be positive. It is this property that allowed Einstein to 
construct a static Universe, balancing the gravitational attraction of normal matter 
with the effective repulsion of the cosmological constant. The total energy density in 
the Universe, Pc can be determined from the expansion rate of the Universe and is 
given by 3H^/{8ttG), where Hq is the Hubble constant. Cosmological data indicates 
that py = 0.7 Pc ~ 4 keV/cm^. A naive estimate of p^ calculated for all modes of 
the zero point field with a high frequency cut-off at the Planck scale is 120 orders of 
magnitude larger than pc- Therefore, there is a need to find a fine-tuned suppression 
mechanism that will bring this large number close to zero. If a small vacuum energy as 
well as a cosmological constant are present, then their values need to be fine tuned to 
the small number p^, which amounts to a fine cancelation to 120 decimal places. This 
bizarre coincidence is the present cosmological constant problem. To find a solution 
to this problem, the contributions from the vacuum fields have to be understood 
and estimated better. Since Casimir force is a direct manifestation of the vacuum 
fluctuations, it provides a tool to comprehend the 'vacuum'. 

1.3.2 "The Hierarchy Problem" 

Interactions in nature have been identified to be of four fundamental type: gravita- 
tional, electromagnetic, strong and weak. General Theory of relativity and Standard 
Model of particles physics are the two most successful theories that explain these in- 
teractions. Attempts to link these two theories are plagued with difficulties due to the 
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vast differences in the strength of gravity as compared to others. The Standard Model 
of electroweak and strong interactions unifies electromagnetic and weak interactions 
at the characteristic energy scale, M^w of 10'^ GeV (= 1 TeV). The energy scale of 
gravity is the Planck energy, Mpi (= 10^^ GeV), where the Compton wavelength of 
a particle becomes equal to its Schwarzchild radius. In natural units {h = c = 1), 
the gravitational constant G = l/M^i. The vast difference between these two en- 
ergy scales is the "hierarchy problem" . Several frameworks have been put forward to 
solve this problem. The most popular among them are string theories, M-theory and 
theories with large extra dimensions (see Plj for a recent review). 

Large extra dimensions and Warped geometries 

A model that attempts to reconcile gravity and quantum theory is one in which the 
fundamental objects that constitute our Universe are not particles but very tiny ex- 
tended objects: strings. However, to have a consistent string theory that can explain 
all known phenomena, spatial dimensions larger than 3 are required (for a review, see 
jl^ HI [I]). The question then is, why do we not see these additional dimensions? 

Two scenarios have been proposed to explain this. One that follows the original idea 
by Kaluza [32], which tries to unify the fundamental forces at Planck scale. In this 
picture, besides the three spatial dimensions of infinite extent, there are additional 
dimensions of finite size, Vc, that are curled up as compactified circular extra di- 
mensions. The typical size of the circle would be determined by the Planck length, 
10^^'^ cm. This means that physics at short distances appears to be higher dimensional 
and forces go as l/r~("+^), where n is the number of extra dimensions. At distances 
much larger than Tc, n = and one observes the law. Around the compatifica- 
tion scale, the forces can be modelled as e''^^'^" and would show up as deviations from 
the standard law. 

The other idea proposed by Arkani-Hamed and others [Sj [3] , tries to unify forces at 
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the TeV scale, thereby removing the hierarchy problem. The 10^^ times weakness of 
gravity as compared to the other forces at this electroweak scale is explained by the 
presence of the extra-dimensions of finite size R^,. They propose that the three spatial 
dimensions in which we live are perhaps just a membrane (3-brane) embedded in a 
higher dimensional bulk of (3 + n) spatial dimensions. The Standard Model fields are 
confined to the 3-brane while gravity can propagate into the bulk. Thus, at distances 
greater than i?*, gravity spreads in all 3+n dimensions and goes as l/r"*^""*"^) while the 
strength of the other forces, still falls as At r > i?*, gravity too reverts back to 

This scenario is distinguished by the term, "large extra dimensions scenario", as 
the additional dimensions in this theory could even be macroscopic. If there were only 
one extra dimension, its size would have to be of the order of 10^° km to account for 
the weakness of gravity. Such an extra dimension would change the dynamics of the 
Solar system and is eliminated by known experimental results. With two equal extra 
dimensions, the scale length would be of the order of a 0.3 mm. This is inconsistent 
with laboratory experiments |22j and astrophysical bounds [2l [121 EOj [21] . For n > 3, 
the scale length is less than about a nanometer. This does not imply that the new 
dimensions will not show observable effects in experiments at sub-millimeter scales. 
A single large dimension of size 1 mm with several much smaller extra dimensions 
is still allowed. Experiments have shown that the scale length of the 'largest' extra 
dimension has to be < 200 fim [221 H]- This would give rise to observable changes 
to inverse square law of gravity at these scales. At such distances, the strengths of 
gravity and Casimir forces are comparable and in order to look for new corrections 
to inverse square law of gravity, Casimir background has to be first understood and 
eliminated. 



New particles 

The super-symmetric extensions to standard model unify the electro-weak and strong 
interactions at energy scales of 10^^ GeV, which is very close to the Planck scale. 
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with the additional assumption that there are no charged particles between the TeV 
and the Planck scale (see for ex. [H |32] and references therein). A trade mark of 
super-string theories is the occurrence of scalar partners of the graviton (dilaton), 
gravitationally coupled massless scalers called moduli, and other light scalers like ax- 
ions. The exchange of these particles could give rise to Yukawa type interactions, that 
would appear in experiments that measure forces. The range A, of these interactions 
depends on the mass of the elementary particle. For super- symmetric theories with 
low energy (few TeV) symmetry breaking, these scalar particles would produce effects 
in the sub-millimeter scales j^SlEO]- The predictions on the strength of these effects 
is less precise than those of the extra-dimension scenarios. 
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Fig. 1.4: Plot of the constraints on the Yukawa interaction parameter a from various 
ranges of A reproduced from Curves 7-10, 12 follow from Casimir force 
measurements. Curve 11 from van der Waals force measurements. Curve 
6 is from an experiment that measured deviations from Newton's law of 
gravity [25j. The typical prediction of extra dimensional physics is shown in 
Curve 13. The region in the (a. A) plane above each curve is excluded and 
the region below each curve is allowed. 



References 



14 



1.3.3 Constraints on new macroscopic forces 

Keeping in mind tliat the strengths of the new macroscopic forces are small, their 
potentials are scaled with respect to the gravitational interaction between two point 
masses as shown below: 

GM1M2 / , , 

V{r) = ^ (1 + ae-^) (1.11) 



where a represents the coupling strength of the interaction and A the range 
The typical scale of A will vary depending on the source of the potential. For the 
extra dimension scenarios it would be the size of the extra dimension, while for the 
new string inspired forces, it would be proportional to the inverse of the mass of the 
mediating particle. Thus constraints can be placed on the parameter space of a — A 
from experiments that study long-range interactions [5]-[ll]. For distances of < 0.1 
mm, Casimir force provides the dominant background and best limits on a for these 
A can be obtained from Casimir force measurements. The available constraints to 
date on a — A from various experiments are summarised in Fig J1.4[ 

The torsion balance experiments described in this thesis are capable of strongly 
constraining theories with macroscopic extra dimensions, apart from measuring the 
Casimir force and its finite temperature corrections. 
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2. TORSION BALANCE- DESIGN AND 
FABRICATION 



Abstract: This chapter is devoted to the description of the main components of the 
experimental set up. An overview of the experimental scheme will be followed by a 
detailed explanation of the various components of the experiment. The procedures 
followed during assembly to reduce systematic and environmental noises will be de- 
scribed. 



2.1 General principle of the apparatus 



Torsion pendulums have been used as transducers for precision measurements for 
over two centuries. They are known for their capabihty to isolate and measure feeble 
effects that would otherwise be difficult if not impossible to observe against the back- 
ground gravitational field of the earth. (see for ex., [Sill]) Our experimental set up is 
aimed at achieving the sensitivity required to measure the finite temperature effect 
in the Casimir force using a torsion pendulum, optical auto-collimator combination. 
[21 H]. The experiment was to be performed in the separation range 2 fim - 10 fim. 
Casimir force at 10 /xm will produce, in static deflection of ~ 10 ^ radians on 

our pendulum. This is well with in the sensitivity of ~ 10~® radians of our optical 
lever. 



The measurement scheme [Fig. 12. Ij broadly consists of a torsional pendulum with a 
flat circular disc, P , as mass element, suspended using a thin strip flbre. The source 
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Fig. 2.1: Figure representing the measurement scheme. L — spherical surface, 
P = pendulum disc, C = coUimating lens, S = source of light, D = Light 
Detector 

of the force field we wish to study is provided by a spherical surface, L , located close 
to one edge of the suspended disc. The torque on the pendulum due to this force 
is inferred from the deflection of the pendulum measured using an auto-coUimating 
optical lever. 

2.2 Torsion Pendulum 
2.2.1 Mass Element 

The mass clement of the torsion pendulum has two parts. A circular disc and a 
mechanism to attach the disc to the suspension. The mass element assembly weighs 
a total of 52 grams and has a moment of inertia of ~ 198 g.cm^ 
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The Disc: 

The active element of the suspension is a circular disc of thickness 4 mm and diameter 
80 mm made of glass. The edges of the disc are chamfered and both the faces are 
polished to have a surface finish of A/2 and are coated with a 1 fim thick layer of 
gold. One of the faces acts as a conducting boundary for the force under study and 
the other face acts as a mirror viewed by a sensitive optical lever which measures the 
angle that the normal to the disc makes with the optic axis of the optical lever. 

The Disc holder: 

The gold coated glass disc is held in a frame made of a gold strip that is 80 fim 
thick and 4 mm wide. In order to hold the glass disc firmly, this strip is shaped 
to form a groove that matches the chamfered edges of the disc and forms a circular 
frame around the disc. The ends of this strip are held securely between the flat 
surfaces of an Aluminium holder [Fig. 12.2] . The flat surfaces that press the strip ends 
together have 50 fim deep, 4mm wide channels machined with the central axis of the 
channel along the axis of the holder. These locate the strip and hence the pendulum 
bob along the axis of the holder. The top of this holder has a 2 mm diameter hole 
through its central axis to hold the torsional fibre. This holder is also gold coated to 
avoid aluminium oxide layers that can accumulate charges. 

2.2.2 The Pendulum Suspension 

The mass element is suspended in two stages to avoid non-torsional modes of oscilla- 
tion of the pendulum. The suspension consists of a pre-suspension, a device to damp 
the simple pendular modes of the fibre, and the main suspension. 
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Fig. 2.2: Picture of Mass element assembly. A = Aluminium holder, P 
Disc 
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Fig. 2.3: A schematic of the pre- suspension mount 
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Pre-suspension: 

The first portion of tfie suspension is a torsionally stiff copper wire of 100 fim diameter. 
Tliis ensures verticality of tlie main suspension that uses a metal ribbon. The pre- 
suspension is essential to avoid spurious torsional effects due to tilts of the suspensions. 
The ends of the stiff wire are passed through a copper ferrule of outer diameter 2 mm 
and bore diameter of about 0.5 mm. The ferrule is crimped such that it holds the 
fibre gently but firmly along its central axis. The length of the Copper wire is 7 cm 
between the ferrules. The wire is mounted from the shaft of a rotary feed-through that 
is attached to the top of the vacuum chamber housing the experiment. To electrically 
isolate the suspension from the chamber, the wire is held through a Macor insulator 
as described below [Fig. 12.3] . The shaft holds a cylindrical clamp made of brass, 
which in turn holds a Macor cylinder. Another cylindrical brass holder is mounted to 
this Macor piece and has a 2 mm diameter hole along its axis. The ferrule at one end 
of the wire is passed through this hole and held in place by a screw. A copper disc 
is suspended from the ferrule at the other end of the wire using a similar mechanism 
[Fig. 12.4] . Kapton insulated copper wires connect the fibre suspension to an electrical 
feed-through mounted on the vacuum chamber. 

Pendular Mode Damper: 

The copper disc (C) passes between the pole pieces of an aluminium-nickel-cobalt 
ring magnets (M) [Fig. 12.4] such that it cuts across the field lines (B) of the magnet. 
The fast pendular oscillations of the fibre and violin modes cause eddy currents to 
flow in the copper disc and dissipate these modes. Thus, these modes will be damped 
out rapidly. Since the copper disc and the magnet are axially symmetric, the very 
slow torsional modes are not damped. 
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(a) (b) 

Fig. 2.4: The pendular mode damper assembly: (a) Shows the general scheme and (b) 
shows a photograph of the assembly 

The Torsional Fibre: 

The main torsional element is a Beryllium-Copper strip fibre of width = 90 fim, thick- 
ness = 9 fim and length = 39 cm. The breaking strength of the fibre is approximately 
100 g. Thus our load of 52 g is well within its breaking strength. Copper ferrules are 
crimped at either end of this fibre as in the case of the torsionally stiff copper wire. 
The fibre is anchored to the copper disc on the top and holds the active mass element 
at the bottom. The fibre is pre-annealed under a load of 50 g at a temperature of 
about 150° C for a day. 

The torsion constant of the suspension kf ^ 0.05 dyne cm rad~^ and the suspension 
has a time period of ~ 406 sec. The thermal amplitude of the suspension in the 
absence of external disturbances is ~ 9 x 10^'' radians at room temperature of ~ 
300° K. 
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Fig. 2.5: The capacitance arrangement for damping torsional modes 

2.3 The Capacitor Plates and Torsion Mode Damping 

The torsional oscillations of the pendulum can be damped to an amplitude of about 
10~^ radians by applying capacitive forces to the suspended disc. Two capacitor plates 
are mounted such that they produce opposing torque on the disc. The net torque 
due to the capacitances is varied by adjusting the individual voltages on them. A 
schematic of the arrangement is shown in Fig. 12.51 A voltage applied to Cl rotates the 
pendulum, P clockwise as seen from above while a voltage on Cc rotates it anticlock- 
wise. By switching the voltages exactly out of phase with the torsional oscillations of 
the pendulum, these oscillations are damped. 

The capacitor plates are made of aluminium and consist of circular plates of 1.2 cm 
diameter that are placed within grounded guard rings. The guard ring is insulated 
from the plate with Macor positioners. The two sets of capacitor plates are mounted 
together on another Aluminium fixture such that they are on diametrically opposite 
edges of the suspended disc and positioned at about 2 mm separation from the disc. 
Shielded Kapton insulated copper leads connect the capacitor plates to electrical feed- 
throughs on the vacuum chamber. The capacitors and the mount are gold coated to 
avoid exposed Aluminium oxide surfaces. 

These capacitors can also be used to perform a null experiment in which the deflection 
of the torsional pendulum is balanced by capacitive forces. The position signal from 
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the optical lever is fed back to control the effective voltage on the capacitor plates 
and the torsion pendulum 'locked' at a fixed position. If the torque on the pendulum 
due to force between the pendulum and the spherical lens surface is modified, voltage 

on the capacitors changes to balance this torque. Thus, the change in voltage on the 
capacitor is a direct signal of the torque acting on the pendulum. 

The varying voltages are generated from a 16 bit DAC in a PCI interface card. 
Typically one volt on the capacitor at 1 mm separation, gives rise to a torque of 
~ 10~^ dyne.cm. The voltage on one capacitor is kept fixed at about 4 V and that on 
the other is varied from V - 10 V so that both positive and negative torques may be 
applied to the pendulum. The information on the angular position of the suspended 
disc obtained from the autocollimator is fed back to a PID loop through software 
(Labview) to control the voltage of the capacitor. When the lens is far away and the 
only force on the suspension is the restoring force from the fibre, the PID loop keeps 
the position of the pendulum locked to about 1/10 of a pixel or 5 x 10~^ radians (1 
sec integration). This is below the thermal amplitude of the pendulum which is of 
the order of 10~® radians. 

2.4 The Spherical Lens and the Compensating Plate 

In our experiment the Casimir force between the suspended disc and the spherical 
surface of a lens is measured. This configuration is simpler to implement as difficulties 
in holding the 2 plates parallel to each other while measuring the force are avoided. 
The lens is 25 mm in diameter and has a radius of curvature of about 38 cm. It 
is coated with 1 /xm thick layer of gold. It is mounted in a cell and held along the 
diameter of the suspended disc close to one edge such that the interactions between 
the lens and the mass element apply a torque on the pendulum. In this position, 
the gravity due to the lens cell assembly will also apply a torque on the suspension. 
This is minimized by using a compensating plate of Aluminium with mass equal to 
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Fig. 2.6: The lens (L), compensating plate (C) and the capacitor plates {Cl and Cc) 
assembled (without the pendulum). T is the translation stage and E the 
EncoderMike® 

that of lens and the cell [Fig. 12. 6j . The outer diameter of the compensating plate is 
equal to that of the lens cell, but the compensating plate has a flat surface facing 
the suspended disc and its thickness is adjusted to equalize the masses. The compen- 
sating plate is mounted such that the gravitational torque due to the plate opposes 
the gravity due to the lens assembly. The net gravitational torque on the pendulum 
is small. More importantly, changes in the force as the lens is moved through small 
distances is negligible compared to the changes in Casimir force and electrostatics 
forces. The lens assembly and the compensating plate are together mounted on a 
translation stage (T) (Newport Model- 461 series) (Fig. 12. 6p . An EncoderMike® ac- 
tuator (E) from Oriel is used to translate this stage perpendicular to the disc surface. 
The actuator movement is controlled by DC voltages applied to it. The encoder has 
a resolution of 0.05 /im and its output is monitored using a Data acquisition card 
with a PCI interface attached to the PC. 

The lens assembly and the compensating plate are electrically isolated from each 
other and from the mount. Kapton insulated copper wires connect them to separate 
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Fig. 2.7: Picture of the torsion pendulum, capacitors and the lens assembled inside 
the vacuum chamber 

connectors on an electrical feed-through attached to the vacuum chamber. 

2.5 The Vacuum Chamber 

The experiment is conducted in high vacuum (~ 4 x 10"*^ Torr). The presence of 
gas surrounding the pendulum at pressures above 10~^ Torr not only damps its os- 
cillations very quickly, but also produces pressure gradients which lead to erratic, 
anharmonic deflections of the pendulum. A cubical vacuum chamber with side ports 
for optical windows with 300 mm long extension tube on its top flange was designed 
and fabricated. The experimental apparatus are mounted inside this chamber [Fig. 12.71 
and Fig. l2.8j . The extension tube is fitted with a rotary feed through at the top, from 
which the pendulum is suspended. This enables us to rotate the pendulum and change 
its equilibrium position. The feed through is motorized and has angular resolution 
of 0.1 degree. The MotorMikes®, translation stages and all other components were 
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vacuum tested to ~ 10 ^ Torr individually before assembly into the main vacuum 
chamber. 

The vacuum chamber is initially pumped down to ~ 10~^ Torr using a turbo molecu- 
lar pump from Varian with a pumping speed of 250 litre/s. An ion pump is switched 
on at this stage. The system is baked at 90° C for 2 days. The turbo pump is then 
valved shut and switched off. The ion pump remains on through out the experiment 
and maintains the pressure at 3 x 10~^ Torr without any mechanical disturbances. 
The pressure inside the vacuum chamber is measured using a compact full range cold 
cathode gauge from Pfeiffer Vacuum. The ion current read out from the ion pump is 
also used to monitor the pressure. 

For the optical lever to view the mirror suspended inside the vacuum chamber a home 
made wire sealed optical glass window is used. The window is made of BK7 glass, 25 
mm thick and is polished on either side to a surface accuracy of A/4. A commercial 
glass view port typically does not have the optical quality required to get sharply 
focused image of the source slits of the optical lever on the CCD detector. There are 
also two other standard glass view ports to see into the chamber. Glass surface being 
dielectric and insulating, accumulates charges. These charges can create electric fields 
inside the vacuum chamber that influences the torsional pendulum. To shield these 
fields grounded wire mesh are mounted inside the vacuum chamber in front of all 
the glass windows. Without these wire meshes, the pendulum can get 'locked' into 
position determined by stray electric fields from the view ports. 

2.6 The Electrical Wiring and Grounding in the Apparatus 

The vacuum chamber is held at ground potential by a thick copper cable. This shields 
the apparatus inside from electrical pick ups. The pendulum and the compensating 
plates are externally connected to the same ground point. The capacitor plates and 




Fig. 2.8: A photograph of the experimental set up during assembly 
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Fig. 2.9: Schematic of the experimental set up 
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the lens are also connected to this ground point when voltages are not applied on them. 

A UV lamp is placed inside the chamber and flashed on during pump down when 
the pressure in the chamber is dropping from 10 Torr to 1 Torr. This generates lots 
of electrons (by photoelectric effect) and some ions; thus allowing a neutralization 
of the electrical charges. The usage of the lamp was 'empirical' and was not very 
systematic. We found small reduction in the residual electrostatic force when the 
lamp was operated for short duration during pump down. 

2.7 The Thermal Panels 

The experimental set up is shielded from fast temperature variations and temperature 
gradients in the environment. The set up is surrounded by four 1.2 x 4.2 meter, 
insulating panels. Each of these panels is made of several layers of thermocol and 
plywood sheets sandwiched between Aluminium sheets and held together by a wooden 
frame. The 'walls' formed by these is covered on the top by thermocol layers attached 
to Aluminium sheet. Various electrical and signal cables come out through tightly 
packed holes in the shroud. The entire apparatus is placed within a closed room 
and controlled from outside. The peak to peak variation in temperature inside the 
enclosure is 1 degree per day (diurnal cycle), while the ambient temperature changes 
by as much as 10 degree per day. However, fluctuations over time scales of an hour 
are within 5 millidegree. 
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3. THE AUTOCOLLIMATING OPTICAL 
LEVER 



Abstract: An optical lever of novel design built to measure the deflection of the tor- 
sional pendulum will be described in this chapter. The optical level has a large dynamic 
range of 10^ and a sensitivity ci/~ 1 x 10"^ radians/^ Hz. The chapter will begin with 
the discussion of the principle of the design and go on to describe its implementation. 
Finally the tests and characterization of the autocollimator will be discussed. 

3.1 Conceptual aspects of the design 

The angular deflection of the torsion balance contains the signal in our experiment. 
The method used for its measurement is the standard optical lever arrangement [see 
for example, [21 E], [1] and references therein] where a beam of light is reflected off 
a mirror on the torsional pendulum and the deflection of the light beam is then 
proportional to the rotation of the pendulum. This basic scheme has been modified 
to give good accuracy and large dynamic range for angle measurements. 

The optical lever is arranged in an auto-collimating configuration. In this configu- 
ration, the translation of the mirror does not change the image on the detector and 
the optical lever is sensitive to only rotations of the mirror. An illuminated array of 
slits, S , is placed in the focal plane of an achromatic lens, C , with a slight offset 
with respect to its optical axis. The coUimated beam emerging from the lens falls on 
the mirror, M , whose rotation angle is to be measured. The reflected beam from the 
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Fig. 3.1: Classical autocollimating arrangement 

mirror returns through lens and will form an image of the slit very close to the slit 
itself, but with an opposite offset with respect to the optical axis. This image falls 
on a linear array CCD detector, D. The location of the centroid of the image on the 
CCD is a measure of the angle between the optic axis and the normal to the mirror. 
In the conventional optical lever, a single slit is imaged on to a photodetector or a 
position sensitive photodiode. The resolution of the device improves roughly as the 
inverse of the square root of intensity of light detected and using a multi-slit scheme 
can increase the sensitivity significantly compared to the single slit scheme. The use 
of the linear CCD array allows measurement of deflections over much larger range of 
angles than what is possible in a conventional optical lever. 

Let i{x) represent the number density of the photons counted as a function of their 
location x. For considerations of the design let i{x) be a simple box shaped distribu- 
tion as shown in Fig. 13.21 
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Fig. 3.2: Distribution i{x) 



The moment, m of the distribution i{x), say about a; = is simply given by 

i{x) ■ {x — Qi)dx; 



m 



X 



(3.1) 



The fluctuation in m is given by, 



/ 



i{x) x^ dx; 



X 



= -(3Z^/; + — ). 



(3.2) 



Thus, the centroid of the image is obtained by dividing the moment by the total 
number of recorded photons Iq = iqw 



m ± Am 
iow 



(3.3) 



Now consider the presence of a background light which generates counts ib{x) spread 
over a width Wb > w about some location l^. This will combine with the image and 
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generate a new centroid given by 

-^0^ + hh I 



Xb 



lo + h 



T ~2 ~T T (3-4) 

{lo + hr +12^+(/o + /,f^'^+ 12^- 

The presence of such a background induces both a systematic uncertainty and an 
additional statistical uncertainty which can be large. 

One element in the design of our optical lever is a strategy to eliminate the error 
due to background light. From the total intensity field it = i{x) + ib{x) we subtract 
+ 3^^^ where p is the width of the digitizing pixel which is much smaller than w 
and Wb- S = 3^^^ is the statistical fluctuation in the background light. After such a 
subtraction we generate a new intensity field given by 

in = {i{x) + ih{x)} -ib + 3i— (3.5) 

V P 

~ i{x) — 36 for in > 0. 

= for in < 0. (3.6) 
The new centroid calculated with this in is given by 



Notice that x„ ~ Xc when the fluctuation in the background intensity, S, are small. 

The second element in the design involves having multiple slits and corresponding 
multiple peaks in the image. Let us consider a grating of u elements with a spatial 
periodicity of w, with ^ opaque and ^ transmitting [Fig. 13. 3j . The full length of the 
grating is W = uw. 

Centroiding as before we get 



3. The Autocollimating Optical Lever 



39 



w 



L-W/2 



L + W/2 X 



Fig. 3.3: Grating 



The uncertainty in Xg is much larger than that in Xc and is about \'2h'Axc- This larger 
uncertainty is essentially due to the fact that the width of the light distribution has 
increased i/-fold. 

To overcome this, consider a set of u fiducial points Xj, i = l,2,...,i/, spaced at 
intervals of w. The centroid of the image of the individual grating elements with 
respect to the corresponding fiducial elements is given by 




(3.9) 



Averaging all the Xc 



(3.10) 



The precision in the determination of the centroid is substantially improved in this 
case as opposed to Eqn. 13. 8[ It is as though the photon density of an individual 
image of the slit has been increased by a factor u. Since, the illumination of the slit is 
limited by the brightness temperature of the source, the above method of decreasing 
the statistical uncertainties proves useful. 



The considerations related to the spacing of the grating are straightforward. The 
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mirror whose deflection angle is to be measured is smaller than the size of the auto- 
coUimator lens and as such determines the diffraction width of the image of the slit. 
Thus, the width of the opaque region between the slits may be taken as twice the 
full width of the diffraction width due to the mirror: 

ao = 2^^/. (3.11) 

Here dm is the diameter of the mirror and / the focal length of the lens. The width 
at of the transparent part of the grating should be chosen such that the width of 
its image covers at least ~ 5 pixels, so that the spatial digitization of the image is 
adequate and does not lead to inadequate sampling of the possible asymmetries in 
the image profile. 

at + ^/>5p. (3.12) 

Although one would like to keep at as small as possible in an attempt to improve 
the angular resolution of the optical lever, it is necessary to keep at to be at least as 
large as the diffraction width itself to achieve a bright enough image, i.e., at ~ 2.5p. 
Further, it may be advantageous to choose its width large enough so that the light 
illuminating it does not get diffracted away even beyond the periphery of the lens. 
Thus, the grating constant a = a^ + at and it is of the order of four to five times the 
diffraction width of the mirror. 

3.2 Construction of the Optical Lever 

A sketch of the optical lever is shown in Fig. 13. 4[ It consists of an array of slits made 
with a photographic plate with transparent slits of 30 /im width, 90 in number, sepa- 
rated from each other with dark regions of 120 /im width (opacity of the dark regions 
is < 100%). This array is illuminated by light emanating from a brightly illuminated 
ground glass sheet placed about five millimeters behind the array. A bank of red 
LEDs which emit in a forward cone angle of ~ 15 degree illuminate the ground glass 
sheet with overlapping circles of light. This increases the brightness and uniformity of 
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Fig. 3.4: Layout of optical lever assembly 

illumination of the slit array. The array is placed in the focal plane of an achromatic 
lens of 1000 mm focal length and diameter 80 mm. The light passing through the 
shts is reflected towards the lens by a right angled prism, located such that the vir- 
tual image of the array is in the focal plane of the lens. The coUimated light passing 
through the lens is reflected by the mirror whose deflection is to be measured. The 
reflected light passes back through the lens and gets imaged on the focal plane, just 
below the virtual image of the source. 

The image is recorded on a CCD camera with a linear array of 6000 Turbosensor® 
photo-elements with a pixel size of 10 /im x 10 /im and center to center spacing of 
10 iim. The well depth is 10^ electrons. Two Analog to Digital Convertors (ADCs) 
built into the CCD camera digitize the charge on each pixels. The even pixels are 
read by one ADC while the odd pixels arc read by the other. A constant bias voltage 
is added to the ADC input to ensure that the ADC always has a non-zero input. 
The camera is controlled using a National Instruments image acquisition card PCI 
1424 and Lab View software. The card sends the clock signals required to transfer 
the charges from the CCD pixels one by one to the ADC and also reads out the 
digitized output of the ADC. The clock signals also control the integration/exposure 
time of the camera. Thus a digitized image of the slit array is obtained. This image 
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Fig. 3.5: A sample image of the slits falling on the CCD 

is diffraction limited by the size of the mirror attached to the balance, yet clear peaks 
corresponding to the slits of the array can be seen, [Fig. 13. 5j . The figure shows the 
counts of the ADC as a function of the CCD pixels. These images were obtained 
with an integration time of 32 ms. The calculation of the centroid from these images 
will be described in §3.31 When the mirror rotates through an angle 6 the image and 
hence its centroid moves through a distance 2F6 on the CCD, where F is the focal 
length of the lens. 



The autocollimator assembly - with the slits, lens and the CCD - is housed such that 
the path from the slits to the lens and back to the CCD can be evacuated down to 
10^^ Torr and stays at about 1 Torr for a couple of days when the pump is off. This 
helps to reduce image wobble on the CCD due to refractive index changes in the air 
path caused by fluctuations of temperature and by air currents. 
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3.3 Implementation of the Centroiding Algorithm 

The centroid of the image is proportional to the angle the normal to the mirror 
subtends with respect to the optic axis. The first step in the determination of the 
centroid is to subtract from the individual pixel values q the dark current and the 
bias. The average value of the sum of dark and bias, di for each pixel is estimated 
from 100 frames without any light falling on the CCD. Accordingly we set the effective 
counts nii as, 

rrii = Ci — di for di < Ci 

rrii = for di > Ci (3.13) 

A typical image profile thus generated would appear as shown in Fig. 13.61 The next 
step in the analysis is to clip off the low intensity background regions from the image. 
The background counts Mb, is estimated as 25% of the average counts per pixel falling 
on the CCD, 

Mb = 0.25 where N = Number of pixels (3.14) 

Mi =mi- Mb for rrii > Mb 

Mi =0 for rrii < Mb (3.15) 

The intensities Mj represent the image of the sequence of the grating slits with the 
background subtracted and the noisy low intensity regions trimmed-off. 



The next step in obtaining the centroid of the distribution takes note of the points 
made in Eqn. 13.81 and Eqn. 13.101 These requires us to establish a sequence of fiducial 
points across the CCD-array so that there is a fiducial point in close proximity to 
each of the images of the grating slits, u in number. The location Xi of the first set 
of fiducials is defined by 

1" 



Xi 



H{i - 1) + 



P 



1,2,3, ...D 



(3.16) 
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where fi is the width of the each sht image in pixels (= 15 in our case) and D is the 
total number of fiducial points. 

Now computing the centroids of the intensity peaks that lie just ahead of each of 
these fiducials and averaging we get the effective location of the image 

L 



X 



S [ Ejipl^ ' 



Notice that in the averaging process we have to include only those values of k for 
which all Mj{j = k to k + fi — 1) are not zero; this is indicated with a prime on the 
outer summation sign. This procedure is similar to folding the image z/-times over 
so that all the peaks line up on each other, x thus obtained is the centroid distance 
modulo fip i.e., the mantissa. This may be added to a constant Xc defined by 

L 

f Xi 



p. 18) 



=1 



to get the complete centroid X: 



X = Xc + x (3.19) 



One final step is needed before we are sure that the best possible centroid has been 
obtained. Consider Fig. 13.61 showing two possible image patterns that may occur for 
the two different orientations of the mirror. 

In order to avoid the enhanced errors that will occur when image profiles straddle the 
fiducial points as in Fig 13.6( b). we need to ensure that the fiducials locate the image 
well within them. To this end, we introduce /i set of fiducials Xj j where j = 1,2 ... fi 
such that Xjj+i = Xij + 1. The image counts at the fiducial locations, Mx^j are 
scanned to locate fiducial sets in which, 

M^^ ^ =0 for i = 1,2,. . .D (3.20) 

The centroid, Xj are calculated as given in Eqns. 13.17113.1^ for all Xj j that satisfy 
Eqn. 13.201 The 'true' centroid, C is chosen to be the one for which Xcj is closest to 
fi/2. 
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Fig. 3.6: Two possible image profiles with respect to the fiducial location, k; profile 
in panel (a) will give an accurate centroid as per Eqn J3.19l profile in panel 
(b) will have larger errors, since the centroid is calculated between k and 
k + n 

3.4 Tests and characteristics of the optical lever 

Several tests were performed to characterize the optical lever. The first was to study 
the dark and bias counts of the CCD pixels. In the CCD the first and last 2 pixels 
are internally shielded from light and are used by the CCD electronics to clamp the 
dark and bias. The bias voltage applied to the ADC convertor, is adjusted such that 
ADC output corresponding to the charges on these pixels is locked at 4-5 counts. 
This keeps a check on drifts due to the internal electronics. Over and above this, 
the counts in each pixel of the CCD were monitored for an exposure time of 32 ms 
without any light falling on the CCD. 100 such points were averaged for each pixel to 
determine the average sum of dark and bias of each pixel. This was later subtracted 
from every frame acquired, to determine the dark and bias subtracted counts per pixel. 

The diffraction limited image of the 30 fim slits falling on the CCD are 5 pixels wide 
and are separated by 15 pixels. For an integration time of 32 ms, the peaks are at 200 
counts compared to the saturation counts of 255 (with 8 bit digitization). The counts 
in the region between the peaks is at 40 as against the dark counts of 5 in the regions 
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of the CCD where no hght is faUing. Several such frames are co-added to increase 
photon counts. Each frame with about 600 counts per peak spread over 5 pixels and 
a total of 90 such peaks leads to a photon shot noise of 2.15 x 10~^ pixels/frame. 
This sets 1.08 x 10^^ radians/frame as the shot noise limit to the angular resolution 
of the optical lever. 

In order to characterize the sensitivity of the system, the optical lever is made to 
look at a rigidly mounted mirror and the centroid of the peak intensity is monitored 
for several hours. The image is recorded every 1.6 sec after summing over 50 frames 
of 32 ms exposure time. A plot of the mean subtracted value of the centroid as a 
function of time is shown in Fig. 13.71 The fluctuations in the value of the centroid, 
as defined by the standard deviation of the data, are found to be at the 3.6 x 10^'' 
radians level, which is worse than the expected shot noise limit of 1.53 x 10~^ radians 
for 50 frames [see Fig. 13.7] . The source of this noise could be the air currents inside 
the autocollimator tube. 




I I I I I I I I I I I I I I I 

1x10' 2x10' 3x10' 4x10' 5x1 0' 6x1 0' 7x1 0' 

Fig. 3.7: Plot of variation in centroid position of a fixed mirror with time. The optical 
lever mount is at atmospheric pressure. Interval between readings is 1.6 sec. 
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To reduce the air currents, the path between the sht and the lens has to be evacuated. 
The fixed mirror is monitored as before after pumping the autocolhmator to 10~^ Torr. 
The variation in the centroid in this case is shown in Fig. 13.81 The fluctuations in the 
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Fig. 3.8: Plot of variation in centroid position of a fixed mirror with time when the 
optical lever mount is at about 1 Torr. Interval between readings is 1.6 sec. 

centroid value are reduced to about 1.28 x 10^^ radians for 50 frames, which is at the 
shot noise limit. Thus we have achieved a sensitivity of ~ 1.6 x 10^® radians/ -^/Hz. To 
improve on this, the total number of photons incident on the CCD will have to be in- 
creased by improving the mask and by efficiently coupling the light source to the slits. 

The dynamic range of the optical lever is ~ 1.8 x 10^. The maximum angle that the 
instrument can measure is defined by the total length of the detector. The present 
mount for the optical lever masks about 500 pixels on either end of the CCD. Thus, 
with a total image width of 1.35 cm, the largest displacement in the centroid that we 
can measure is 3.65 cm, which corresponds to 1.8 x 10^^ radian defiection of the mirror. 

Even at this stage we can see the salient features of this design of autocollimating op- 
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tical lever: it is capable of measuring angles with a sensitivity of ~ 10~^ radians/v'Hz 
and has a dynamic range exceeding a few millions. These features are quite unique. 
With a better mask and hght source we hope to reach a sensitivity of 3 x 10~^ 
radian/ ^Hz without sacrificing the dynamic range. 
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4. CHARACTERIZATION OF THE 
APPARATUS 



Abstract: This chapter presents the tests and characterizations performed on the ex- 
perimental apparatus. 

4.1 The Torsional Pendulum 

The force transducer in our experiment is a torsional pendulum. The sensitivity of 
the torsional pendulum to any torque acting on it is related to the moment of inertia, 
/ of the mass element suspended and torsional constant of the suspension fibre, n. 
These parameters determine the time period of oscillation of the pendulum. The 
period, T;, is given by 



Another important parameter of an oscillator is its quality factor, Q. It is a measure 
of the damping in the system and is defined as the ratio of time taken for the ampli- 
tude of the oscillator to fall by 1/e to the time period of the oscillator. 

In order to characterize our torsional pendulum, it was allowed to oscillate freely inside 
the vacuum chamber at a pressure of 4.0 x 10~^ Torr with an amplitude of ~ 3.0 x 10~^ 
radians. The oscillations are monitored continuously for about a day. Fig. 14.11 shows 
a few cycles of the oscillations. The time period and Q are obtained from by Fourier 
power spectrum of this time series data. The power spectrum calculated from the Fast 
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Fig. 4.1: Free Oscillations of the pendulum 
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Fig. 4.2: Power spectrum of the torsional oscillations 
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Fourier Transform(FFT) of the data is shown in Fig. I4.2[ The peak in the spectrum 
occurs at a frequency of 2.44 x 10"^ Hz corresponding to a period of 409 sec. The 
quahty factor as determined from the width of this peak is ~ 56. But this is hmited 
by the resolution of the FFT fixed by the total observing time of about 20 hrs. The 
true Q is ~ 100 as we can see from several such one day stretches of oscillations. The 
equilibrium position of the pendulum drifts by ~ 5.0 x 10~^ radians, peak to peak, 
over a day (diurnal cycle). The signal to noise ratio achievable in the system depends 
on the Nyquist noise which puts a limit on the smallest torque that we can measure. 
The Nyquist noise is given by. 



Here, k is the torsion constant of the suspension fibre, T the laboratory temperature, 
Q is the quality factor of the fibre, t is the total time of observation and ujh is natural 
oscillation frequency of the balance. Taking typical values T = 300° K, k = 0.05 
dyne.cm.rad~\ Q = 100, t = 50 s, and = 1.5 x 10~^ rad.s"\ r^yg ^ 6.1 x 10~^ 
dyne. cm. This value is two orders of magnitude smaller than the force sought for in 
the present set of experiments. 

4.2 Pendulum behaviour during evacuation 

It is repeatedly noticed that the level of vacuum in the chamber plays a crucial role 
in determining the trajectory of the pendulum within the assembly. During assembly 
the lens is typically positioned at a separation > 5 mm from the pendulum while the 
capacitors are fixed at about 2 mm from it. When pump down starts, the pendulum 
tends to remain 'stuck' to either of the capacitors until a pressure of ~ 10^'^ Torr 
is reached. At this stage the pendulum starts to oscillate with very short periods 
of about 20 s. As the pressure drops the period progressively increases and reaches 
~ 200 s at 10~^ Torr. The period stabilizes to the 'free oscillation period' of about 
400 s when the pressure reaches ~ 10"'' Torr. Further drop in the pressure does 
not change the period of the pendulum. We have also tested the behaviour of the 




(4.1) 
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pendulum when the pressure is allowed to rise due to degassing, with pumps shut 
off. The period begins to drop as the pressure increases to about 10~^ Torr. This 
behaviour of the pendulum was not systematically characterized. But, it seems clear 
that even the low pressure cushions sandwiching the pendulum disc in the small gaps 
between the capacitors and the disc can severely alter the effective torsional spring 
constant of the oscillator. Our system is operated at a vacuum of ~ 4 x 10~® Torr to 
avoid these effects. 



4.3 Calibration of the force arm 



The experiment attempts to measure the Casimir force between a spherical lens sur- 
face and the fiat surface of the torsion pendulum. The torque due to this force on 
the pendulum, depends on the lateral distance (force arm) of the point at which the 
force acts on the pendulum as measured from the suspension axis. The force arm 
is determined from the point of contact between the lens and the pendulum, which 
is estimated as follows: at some initial position of the lens, the pendulum is pulled 
towards the lens by applying voltages on the capacitor plates. The pendulum position 
is monitored using the optical lever as it approaches the lens. The position where 
the pendulum touches the lens is deciphered from the centroid location on the optical 
lever at which the pendulum suddenly turns around. The average point of contact 
of the lens is measured from several touches at each position of the lens. The lens is 
then translated by known amount using the EncoderMike® actuator [see §2.4]. The 
point of contact of the lens on the pendulum and the force arm change negligibly, 
but the angular position of the pendulum at which the contact occurs changes signif- 
icantly [Fig. 14. 3j . This new angular position at which the contact occurs is measured 
for several displacements of the lens and a plot of it is shown in Fig. 14. 4[ Thus, the 
relation between the linear displacement of the lens and the angular displacement of 
the centroid of the image at the point of contact is obtained from the slope of the 
data. The slope tells us that 0.114 fim translation of the lens corresponds to 1 pixel 
movement of the image centroid. Hence when the point of contact moves transverse 
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Fig. 4.3: Schematic representing the 'touch position' and its relation to the force arm 
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Fig. 4.4: Plot of the touch location as a function of lens position 
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to the plane of the pendulum (towards or away from the lens) by 0.114 fim, image 
shifts by 1 pixel. From the optical lever parameters, we know that a image displace- 
ment of 1 pixel corresponds to a transverse movement of 0.2 fim at the edge of the 
pendulum disc which is at d/2 = 40 mm from the suspension axis [Fig. 14.3] . Thus 
the force arm / or point of contact is at 22.8 mm from the axis of the pendulum. 



5. CASIMIR FORCE MEASUREMENT - 
STRATEGY, DATA ACQUISITION AND 
ANALYSIS 

Abstract: The basic procedures followed in our experiment will be described in this 
chapter. The chapter begins with an overview of the forces influencing our torsion 
pendulum transducer and goes on to describe how the data is acquired in order to 
comprehend these. The analysis of this data to extract information about Casimir 
force is then discussed in detail. There is clear evidence for the finite temperature 
Casimir force, detected for the first time, in the analyzed data. 

5.1 Strategy of the experiment 

The main objective of the experiment is to measure the force of attraction between a 
fixed lens and a fiat plate suspended as a torsional pendulum. In the absence of the 
lens, the only torque acting on the pendulum is the restoring torque from the fibre. 
As the lens is brought close to the fibre, various forces begin to act on the pendulum. 
The basic torques acting on the pendulum when the pendulum oscillates near the 
lens are the fibre restoring torque and the torques due to electrostatic and Casimir 
force. 

The presence of all these forces will introduce changes on the torsional oscillations of 
the pendulum. Since these forces in general arise from potentials that differ from the 
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quadratic harmonic potential due to the torsion in the fibre, they make the motion 
of the pendulum anharmonic. As described in the previous chapters, we record the 
angular position, 9 of the fibre as a function of time, t. Thus the strategy adopted is 
to extract information about the forces from the t versus 9 data. 



The free oscillations of a pendulum will be simple harmonic. The existence of other 
potentials, besides the harmonic potential due to the fibre, will contribute to the 
acceleration 9, thereby changing the simple linear relationship of 9 with the amplitude 
of the oscillations. If f{9) is the net force on the pendulum then. 

The torque on the pendulum, 

t{9) = f{9)l and (5.1) 
l'9 = t{9) (5.2) 

where I is the lateral distance of the point at which the force acts on the pendulum 
as measured from the suspension axis, 9 is the angular acceleration of the pendulum. 
Similarly the angular velocity 9 is related to the potential U{9), 

h9^ + U{9) = const. (5.3) 

where t(9) = (5.4) 
o9 

Thus, by measuring 9 or 9 the potential and the forces that affect the pendulum can 
be characterized. 



5.1.1 Summary of Forces acting on the Pendulum 

Several forces between the lens and the pendulum begin to become sizeable as they 
approach each other to within about 100 //m. We now discuss the forces that play a 
major role in the motion of the pendulum. 
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Fibre restoring torque: 

The restoring torque due to the fibre makes the pendulum execute simple harmonic 
motion about a mean position, hereafter referred to as the 'null' position Oq. As is well 
known, this torque always acts towards 60 and in proportion to the defiection of the 
torsion balance from this position. In the following descriptions, the distances and 
the angles are measured with respect to the position of the lens as origin [Fig. 15. Ij . 
During the experiment, the presence of other forces will generate new torques which 
will shift the balance to a new equilibrium position O^.. By appropriate rotation of the 
shaft from which the balance is suspended, the equilibrium position 6^ is made to lie 
between the lens at 9i = 0, and Oq. The restoring torque, Tfn, and the acceleration, 
9fih on the pendulum due to this torque are given by: 

= f^iOo-e) (5.5) 

9fib = j{9o-0) (5.6) 
Ofib = 2.46 X 10-^(^0 - 0) (5.7) 

where k is the torsion constant of the fibre, I the moment of inertia of the pendulum 
and 6*0 is the equilibrium position of the pendulum. Similarly the potential, Ujib and 
the angular velocity, Ofn, are given by, 

Uf,t = '^{00 - ey + const. (5.8) 

It is convenient to write in general, 

2U 

9ji^ = — — + const. (5.9) 
= -Un + const. (5.10) 

With this we may think of 6"^ as just the negative of the normalized potential 11^. 
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Fig. 5.1: Reference axes for the system. 



Casimir Force 



The Casimir force per unit area between two parallel conducting plates separated by 
a distance at T = K is given by, 



fco{z) 



240^4 

0.013 _ . . 
— dyn.cm where = z m microns 



(5.11) 
(5.12) 



For a finite temperature T we approximate the distance dependence as 



fcTiz) 



fco{z) 



for z < \t 



fcxiz) = /co(2)t— for z> Xt 
At 



(5.13) 
(5.14) 



where At is the distance at which the finite temperature effects become important. 
Comparing Eqn. 15.141 with the expression for Casimir force at 'high temperature' 
described in Chapter 1 (Eqn. 1.5), 



Fj(rf)^-i:^^ at high T 



(5.15) 



we can estimate that for T = 303° K, ~ 3.26 /xm. 
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L[nm] 



Fig. 5.2: Ratio of the Casimir force calculated taking into account the finite conduc- 
tivity and finite temperature effects to the Casimir force at zero temperature 
between ideal conductors is plotted as a function of spacing between plates 
(from, [IT]). Our experiment is in the distance range of 2 fim - 10 /im 

For convenience, we write Eqns. 15.111 and 15.141 as, 

fco{z) = J dyn.cm''^ (5.16) 

z 

fcAz) = dyn.cm''^ (5.17) 

The Casimir force between the lens and the pendulum can now be estimated using 
the "Derjaguin (proximity) approximation" , often called "proximity theorem" p!]. For 
our geometry the expression has two parts to it depending on the separation between 
the lens and the plate. At small separations, the zero temperature Casimir force 
dominates. As the separation increases, finite temperature effects become important 
and the force law changes. Fig. 15.21 shows this change in the law (see chapter 1 for the 
definition of the ratios). The Casimir force for the lens-plate configuration is given 



5. Casimir Force Measurement - Strategy, Data Acquisition and Analysis 



60 



by (see Appendix B for details), 
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for z > Xt 
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where Xmax = ^^ff^ and z = 19 is the separation between the lens and the plate, 
R is the radius of curvature of the lens and r^ax is the aperture radius of the lens 
(Fig. 15. ip . The potentials for the same configuration are. 
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(5.19) 



(5.20) 



3Xt \z {z + Xmax 

The angular velocity and acceleration due to Casimir force can be calculated using 
Eqns. OandEai 



It may be noted that the roughness of the surface is not a major correction term to the 
Casimir force in our experiment, since it focuses on Casimir force at relatively large 
separations between the metal surfaces. However, the finite conductivity correction, 
a reduction in force amounting to approximately 10% needs to be applied throughout 
the measurement range. 

Electrostatic Forces 

There are various sources of electrostatic force between the lens and the plate. Several 
effects lead to potential differences between metal coated surfaces even when they are 
externally connected to ground [121 [131 IIJ • Small difference in the contact potentials 
in these connections can lead to potential differences between these surfaces. With 
care, such potentials can be minimized to a few mV level. Another source of electro- 
static interaction between the bodies are patch fields. These are microscopic electric 



5. Casimir Force Measurement - Strategy, Data Acquisition and Analysis 



61 



potentials found on surfaces of metals. The different crystalline planes of a given ma- 
terial can have work functions that can differ by as much as 1 V . If the coated metal 
surface is a mosaic of random microscopic crystal planes, local potential differences 
will occur between each of these micro-size patches. Surface contamination in the 
metal coatings can also lead to local variations in the work function and also provide 
sites for trapping electrical charges. 

Force due to Capacitative coupling: 

The conductive surfaces of the lens and the plate, form a capacitance which is given 
by "proximity theorem" to be. 



where eo is the dielectric constant in free space, R is the radius of curvature of the 
lens, z the separation between the lens and the pendulum in metres and Vmax is the 
aperture radius of the lens. If the voltage difference, V between the lens and the plate 
were zero, the net force due to this capacitance will be zero. But potentials can be 
present due to any of the reason mentioned above. The electrostatic potential energy 
due to these between the lens and the plate is given by. 



Ci = 




(5.23) 



(5.21) 



(5.22) 




(5.24) 



(5.25) 



The resulting capacitative force is given by. 



nReoV^ 




(5.26) 



(5.27) 



z 



(5.28) 
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We approximate the corrections to the above formula, due to surface roughness and 
a variety of other causes and write for the torque, 

(5.30) 

where the roughness parameter 6r will be determined self-consistently from the data. 
The voltage V comprises of the applied voltage Vap, the D.C voltage V(o,dc) due to 
contact potentials and other causes and V(o,ac) due to the stray pick up fields and 
leakage, etc. This contributes to the normalized force a term. 

Force due to charges on the lens/pendulum: 

Stray electric charges on the surfaces of the lens and the plate also lead to electrostatic 
forces. The potential between these charges depends in detail on their distribution 
on the surfaces and cannot be modeled easily [see for example, [IH]]- An upper 
limit to the strength of this force may be obtained by assuming that all the charge 
is concentrated at the point of contact between the lens and the plate and that a 
simple inverse square law applies for the force with an equidistance image charge. 
The potential due this is given by, 

U,, = (5.32) 
Zm = 2z, the separation between the lens and the pendulum. (5.33) 

The force and the acceleration on the pendulum due to it are. 



Torque, r„ = -^g^, (5.35) 
= -2.94 X 10^^^. (5.37) 
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where q is effective charge on the lens/pendulum. For a distributed charge the ~ 1/0^ 
dependence will change over to a slower dependence of ~ 1/9. This will be the case 
with metallic surfaces when there are some residual charges on them. 



Force due to charge and voltage on the lens: 

One more source of electrostatic force is the combined effect of the charges and voltage 
differences present. The electric field, due to these voltages, V interacts with the 
charge, q to apply an additional force on the pendulum of the form, 

A, = -q£ = -^, (5.38) 

z 

Torque, r^g = —- — , (5.39) 
qV 

Acceleration, 9 = —-fjr, (5.40) 

19 

= -5.2 X 10^4-- (5-41) 
9 

= V(o,dc) + V{o,ac)', however, since the acceleration is linear in V as shown in Eqn 15.411 
and the pick up voltages are at frequencies much higher than the frequency of oscilla- 
tion of the balance at a few milliHertz, the A.C torques will average to zero and will 
not contribute. 



Gravitational force: 

The gravitational attraction of the lens on the pendulum gives rise to a torque of 
~ 4 X 10~^ dyne. cm on it. This corresponds to the attraction between disc shaped 
masses of approximately 10 g each separated by an effective distance of about 5 mm. 
The 'compensating plate' (refer §2.5) is used to balance this gravitational torque 
of the lens. The mass of the compensating plate is equal to that of the lens but 
its position with respect to the torsion pendulum disc is not exactly the same as 
the lens. Therefore, the cancellation of the gravitational torque is only approximate 
and the residual gravitational torque on the pendulum is estimated to be ~ 10~^ 
dyne. cm. This changes only by ~ 10^® dyne. cm when the lens is moved over the 
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distance of 50 /xm, since the fractional change in the effective separation between the 
centers of gravity of the plate and the lens assembly is only 1/100. The change in 
angular acceleration of the pendulum due to this is below the present sensitivity of 
the experiment. 



7 F F F T F 1 F r ' '' ' r ' ' ' r 




# 0.002 0,004 0.006 0.008 

9 (radians) 

Fig. 5.3: A plot of normalized potential f/jv = 2U/I is shown as a function of 6. 

The solid line maps the potential due to the fibre restoring force when the 
equilibrium position of the pendulum, 60 = 5.3 x 10"^ radians. The dotted 
line is a plot of the net potential due to the lens and the fibre for V = 5 mV, 
g = 1 X 10^^^ C and A = 3.26 fim. The dashed line is a plot for V = 50 
mV and the dashed dot line for V = 75 mV on the lens 

5.1.2 The Effect of these Forces on the Pendulum 

A plot of the square of the angular velocity, 6'^ = —U^ of the pendulum as a function 
of 9 is shown in Fig. 15. 3[ This maps the potential in the region of oscillation of the 
pendulum in the absence of the lens (solid line) and when the lens is at a distance 
of about 120 /xm from the null position 6q (various dashed lines). For calculating 
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the potential due to the lens, the net voltage difference between the lens and the 
pendulum V is set to 5 mV (dotted line), 50 mV (dashed line) and 75 mV (dashed 
dot line). The net charge, q is taken to be 1 x 10~^^ C and the distance at which the 
Casimir force changes from zero to finite temperature form, A is taken to be 3.26 /xm. 

Prom the plot it can be seen that the presence of the electrostatic and Casimir forces, 
considerably changes the potential between the lens and the plate. The potential is 
highly attractive towards the lens for small angular separations. The angular distance 
at which this change over to attractive potential happens depends on the net voltage 
on the lens. Thus, when the voltage on the lens is 5 mV, the change over happens 
at a distances of about 5 x 10~^ radians. As the voltage on the lens increases, the 
change over point extends away from the lens. The equilibrium position of oscillation 
of the pendulum is also shifted for voltages higher than 20 mV. The potential beyond 
7 X 10~^ radians is shifted for higher voltages but is parallel to the pontetial due to 
free fibre. The potential below 7 x 10"^ radians is flatter than the potential due to 
free fibre. Thus the oscillation period of the pendulum is modified even for small 
amplitude oscillations of the pendulum. 

As the lens is brought close to the pendulum the signatures of the electrostatic forces 
can be seen by looking at changes in: 

(a ) the equilibrium position of the pendulum, 

(b) the time period of the pendulum, and 

(c) the angular acceleration 6 of the pendulum as a function of its angular position 6. 

Since the Casimir force falls off more rapidly than the electrostatic force, to study 
the effects of Casimir force, we need to go closer to the lens, ie., < 1 x 10~^ radians. 
At those angular distances, the pendulum will no longer oscillate but will start to fall 
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towards the lens. The force in that region then can be elucidated only from the 6 
versus 6 of the pendulum. Hence, to estimate Casimir force, we need to look at the 6 
versus of the pendulum at angular distances of < 1 x 10~^ radians. 

To completely fathom the electrostatic forces though, in addition to measurements of 
the changes on the pendulum as a function of lens position, knowledge about these 
changes as a function of voltage on the lens are also essential. The force due to the 
capacitive coupling increases as the square of the net voltage on the lens, while the 
charge effect is linear in voltage. Hence, by studying the force dependance on the 
voltage, an estimate of the effect of both the residual voltage and the charge on the 
pendulum can be obtained. 



5.2 Data Acquisition: 

Despite considerable care taken during the assembly to minimize the voltages and 
charges on the lens and the pendulum, certain residuals are unavoidably present in 
the apparatus and it is essential to characterize these electrostatic forces. All the 
earlier experiments that looked for Casimir forces have seen the effect of these forces 
[HI m [ISl [IHl [211 [in [H [2]- Since the distance range in our experiment is larger than 
the earlier experiments on Casimir force, electrostatic forces are a strong background. 
Hence, the data acquisition was tuned to be able to estimate the electrostatic forces 
and subtract them in order to look at the Casimir force at separations between 2 fim 
- 10 fim. 

The Casimir force falls off more rapidly with distance than the electrostatic forces. 
Thus, at distances > 25 fim the electrostatics forces will dominate by more than two 
orders of magnitude. At these separations corresponding to a angular distance of 
> 1 X 10~^ radians, even for a net voltage, V of 75 mV between the lens and the 
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pendulum, it will execute anharmonic oscillations. Hence, the presence of the elec- 
trostatic force can be studied by looking at the oscillations of the pendulum in the 
potentials summarized in the previous section. 

Data was acquired in three main modes: 

• Small amplitude oscillations that sampled the region within 5 x 10"'^ radians of 
the equilibrium position. 

• Large amplitude oscillations that sampled regions up to an angular separation 
of 1 X 10^'^ radians from the lens. 

• Accelerated fall of the pendulum onto the lens starting from an angular sepa- 
ration of < 5 X 10~'^ radians. 

Small amplitude oscillations 

To start with, the lens is positioned at a distance of about 0.35 mm from the equi- 
librium position of the pendulum. The oscillations of the pendulum for amplitude of 
~ 2.5 X 10^'* radians, i.e, 5.7 /im are monitored using the autocollimator. The lens is 
then moved stepwise slowly towards the pendulum until the separation between the 
two is reduced to 100 /im. At each separation, the pendulum is allowed to oscillate 
with an amplitude of ~ 2.5 x 10~^ radians for several (> 10) cycles. The position 
of the pendulum in units of CCD pixel numbers is recorded every 0.16 sec. Fig. 15.41 
shows the raw data of few cycles of oscillations of the pendulum. 

From this data the change in the period and the equilibrium position of the pendulum 
as a function of lens position is derived. As the lens is brought closer the equilibrium 
position, the period of the pendulum remained unchanged until the separation is 
> 250 fim. For lower separations, the equilibrium position of the pendulum shifts 
closer to the lens. The period changes from about 406 sec at 0.35 mm separation to 
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Fig. 5.4: The raw data showing a few cycles of small amplitude oscillations of the 
pendulum 



about 503 sec at 100 /im separation. 



These changes are shown in Figs. 15.51 and 15. 6[ Fig. 15.51 shows the change in equi- 
librium position as the lens is moved closer to the pendulum. The pixel numbers on 
the CCD are used as the reference to see this change. The equilibrium position is not 
measured with the lens as origin as the position of lens itself changes. The separation 
marked is the effective separation calculated taking into account both, the change in 
equilibrium position and the change in the lens position. Fig. 15.61 shows the variation 
in the time period of the pendulum for the same separations. 

Thus, at separations > 250 the interaction of the pendulum with the lens is neg- 
ligible and the motion of the pendulum is defined by properties of the suspension 
fibre. Hence, this region can be used to characterize the fibre properties. To study the 
effects of the interaction with the lens, separations < 120 fim are ideal as the change 
in the equilibrium position and time period are steeper in this region. 
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Fig. 5.5: The change in the equihbrium position of the pendulum as the lens is moved 
closer to it. 
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Fig. 5.6: The change in the time period of the pendulum as the lens is moved closer 
to it. 
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In addition to looking at the changes in the equihbrium position of the pendulum 
as a function of separation, the changes in the equilibrium position are also studied 
as a function of voltages applied on to the lens. The lens is kept at a fixed position 
and the voltage on the lens is scanned until the pendulum sees the least force from 
the lens. This would be the voltage at which the applied voltage Vap balances the 
residual voltage on the lens Vo.dc- The data is taken for different null positions of the 
lens. The observations are summarized in Fig. 15.71 
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Fig. 5.7: The equilibrium position of the pendulum as a function of voltage applied 
on the lens at various lens positions 



In Fig. 15. 7[ the separations marked are the distance between the lens and the pen- 
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dulum equilibrium when the apphed voltage on the lens is zero. The separation of 
117.59 /im corresponds to the data taken when the lens is at pixel number 1417, 
104.77 iim to lens at pixel number 1500 and 88.5 /xm to lens at pixel number 1585 
respectively. As the voltage on the lens is increased, the equilibrium position shifts 
away from the lens, implying a reduction in the force of attraction from the lens until 
the voltage reaches between 55 to 60 mV beyond this voltage the attraction increases 
again. Hence, the electrostatic forces are minimum around this voltage range. 

Thus, from the small amplitude data, we have an idea of the magnitude of the elec- 
trostatic forces acting on the pendulum. The distance and the voltage ranges over 
which observations are to be carried out in order to fully characterize the electrostatic 
and other forces are now clear. 

Large amplitude oscillations 

Another mode in which data is obtained is by sampling the potential in the regions 
that are very close to the point where the potential becomes completely attractive. 
As we observed in the small amplitude data, the electrostatic force on the pendulum 
is the least when the voltage is about 57 mV on the lens. This then means that at 
angular distance of about 1 x 10~^ radians, the potential in the region between the 
lens and the pendulum would be considerably different from the harmonic potential 
of the fibre, but will still have a minimum. Thus, if the pendulum is made to exe- 
cute oscillations of very large amplitude about 2 x 10~^ radians, the oscillations of 
the pendulum will show marked anharmonicity especially during one half of the cycle. 

The period and the equilibrium position of the pendulum in this case will be different 
from those in the small amplitude case. The 9 and 9 of the pendulum will show marked 
change between the two halves of the cycle. Since in this mode we are sampling the 
potentials in regions that are reasonably far away from the lens - separations > 20 /im, 
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Fig. 5.8: The raw data showing large amphtude oscillations of the pendulum when an 
voltage of 75 mV was applied to the lens 

the dominant force on the pendulum will be the electrostatic force. Thus, quantitative 
information about the electrostatic forces can be obtained from this data. 

The amplitude of the pendulum is increased using the capacitor plates that were 
described in Chapter 2. The amplitude of the torsional oscillations of the pendulum 
can be controlled by applying appropriate voltage on these plates. The pendulum is 
made to oscillate such that it approaches to within 1 x 10^'^ radians from the lens for 
about 10 cycles. Data is taken at 0.16 sec time intervals. The time period of these 
oscillations varies from 400 sec to 600 sec depending on the voltage present on the 
lens. The angular position of the pendulum is recorded in units of pixel numbers of 
the CCD. Fig. 15.81 displays a portion of the large amplitude data taken and shows 
the asymmetry in the oscillations on the side closer to the lens (lens is at about pixel 
number 1580 for the data shown). To fully understand the electrostatic forces, the 
large amplitude data is obtained by applying voltages on the lens. Since we expect 
the minimum in the force to be at about Vap = 57 mV, voltages starting from 50 mV 
to 100 mV are applied on the lens in 5 mV steps. 

Further analysis of these oscillations to comprehend the electrostatic forces will be 
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presented later in this chapter. 
Pendulum 'falling' onto the lens 

The Casimir force falls very rapidly with the distance from the lens. In order to 
measure Casimir force, the potential very close to the lens has to be studied. But in 
this region, the potential is highly attractive and the pendulum cannot oscillate. In- 
stead, it is pulled steadily towards the lens until it hits the lens and bounces away from 
it. Thus, to quantify Casimir force, data is taken as the pendulum falls on to the lens. 

The initial velocity of the fall is controlled by using the capacitor plates that damp 
the angular oscillations. By applying appropriate voltages to the capacitor plates, 
the pendulum is first pulled over the potential hill onto the attractive region of the 
potential. As the pendulum falls towards the lens, it is stopped and mildly pulled 
away from the lens. The velocity away from the lens is not sufficient for the pendulum 
to go over the potential hill again, but it just falls back onto the lens after moving 
away from it for a short distance. Thus, the fall data is available with zero initial 
velocity. [Fig J5.9j 

Data is taken by making the pendulum fall onto the lens several times. The angular 
position of the pendulum 6 is recorded every 0.16 sec as the pendulum falls. The 
voltage on the lens is also varied as in the case of the large amplitude oscillations. 
Fig. 15.91 shows the raw data of a single 'fall' of the lens on to the pendulum when 
a voltage of 80 mV is applied to the lens. The analysis of this data to extract 
information about the Casimir force is presented in the next section. 

5.3 Analysis: 

The angular position data obtained from the autocollimator, has to be analyzed to 
gain knowledge about the forces acting on the pendulum. The strategy followed 
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Fig. 5.9: The raw data showing the pendulum 'faUing' onto the lens after being pulled 
away from it using the capacitor plates. For this data the lens is at pixel 
number of about 1580 and a voltage of 80 mV is applied to it. 



involves the use of ^ as a function of 9 to extract information about the forces. Thus, 
the first step in the analysis is to convert the {t — 9) data to {9 — 9) map. This is then 
further analyzed to extract information about the force. The analysis of the data 
acquired as described in the previous sections is presented in this section. 



Large amplitude oscillation: 

The large amplitude oscillations of the pendulum give us information about the forces 
acting on the pendulum at separations > 20 yum from the lens. The electrostatic forces 
are the dominant forces at these separations. 



An angular position plot of large amplitude oscillation is shown in Fig. 15.81 As a 
first step in the analysis of the data, the angular position is converted to angular 
separation, 9, taking into account the lens position during the run of the experiment. 
This data is filtered to remove 9 values that have > Sa errors in them. These errors 
arise due to the non-uniformities in the CCD pixels. The 9 values with errors > 3(T 
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Fig. 5.10: Plot of smoothened t — 9 Data taken with a voltage of 75 mV on the lens 

are replaced by the mean 9 of the neighbouring 10 points. Since the time period of 
the pendulum is much larger than the sampling time between points, the change in 
9 within this 10 point interval of 1.6 sec can be considered linear. 

The {t — 9) data thus generated is further smoothened using a third order polynomial 
fit for every 51 points. This is a small fraction of the total of about 2500 points 
available per cycle of the pendulum. The central, 26*'^ point of this set is replaced 
by the smoothened 9 obtained from the fit. A plot of the smoothened data is shown 
in Fig. 15.101 The next step is to obtain the angular velocities and accelerations. For 
this, in every 51 points of the smoothened {t — 9) data, the time at the 26*^^ point ^26? 
is subtracted from all te. This is then fitted with a second order polynomial in t. Let 
this he u + vt + wt^. Here u would be 6' at t = which is the 9 of the 26*^^ point. The 
angular acceleration and velocity are then calculated as follows. 



9{t) 



u + vt + wt^, 



(5.42) 



9 



V + 2wt, at 9 = u, 



(5.43) 



9 



2f , at 9 = u. 



(5.44) 



This is done for entire {t — 9) data by sliding over every consecutive 51 points. Thus, 
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Fig. 5.11: {6 — 6) for large amplitude oscillations data taken with a voltage of 75 mV 
on the lens and a fit to the data 



we get 6, 9 at every 9. Since the (t — 9) is cyclic, the {9 — 9) and {9 — 9) data will also 
be cyclic. The cycles are overlapped by ordering according to increasing 9. These are 
then averaged over a 9 range of 1 x 10~^ radians. Thus, for every set of {t — 9) data, 
a set of 9, 9, 9 data is generated. 



Fig. 15.111 shows a plot of {9 — 9) data. The 'asterisks' mark the data points obtained 
as described above. Almost periodic errors are present even after the averaging as 
can be seen. The error is large when the velocity of the pendulum is large and is 
perhaps caused by systematic errors in the timing generated by the PC. The PC runs 
Windows operating system and hence interrupts from the operating system could 
lead to non-uniformity in the time difference between 2 successive data acquisitions. 



The linear portion of the plot, for 9 > Ax 10^'^ radians represents the force due to the 
fibre. The sharp curvature in the plot for 9 < 2x 10~^ radians is predominantly caused 
by the electrostatic forces. To study these forces, described in §1 of this chapter, we 
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fit an equation of the form, 

The coefficient h represents the force due to the fibre and c the voltage dependant 
electrostatic forces. The coefficient d is an estimate of the Casimir force at these 
separations. It was seen that the fitting coefficients Or and for this data at large 
separations, do not play a major role in determining the best fit parameters. The 
values of coefficients a, h and c do not change much even if 9^ and d are made zero. 
Parameters a and c change by about 0.001% while parameter h does not change at 
all. A fit of this form to the [9 — 9) data when a voltage of 75 mV is applied to the 
lens is shown in Fig. 15.111 



Tab. 5.1: Table of coefficients of the fit to — for large amplitude oscillation data 
obtained with various voltages applied to the lens 



Voltage on lens 
(volts) 


Fit parameters 


a 


b 


c 




d 


0.05 


1.278 X 10-6 


-2.3 X lO-'^ 


-5.0136 X 10-1" 


-6.3578 X 10-^ 


-1.198 X 10-" 


0.055 


1.2433 X 10-6 


-2.3 X lO-'^ 


-4.2532 X 10-1" 


-6.7919 X 10-* 


-1.5841 X 10-13 


0.06 


1.2636 X 10-6 


-2.3 X lO-'^ 


-4.4659 X 10-1" 


1.3246 X 10-* 


-9.2951 X lO-i-i 


0.065 


1.2612 X 10-6 


-2.3 X lO-'^ 


-4.0365 X 10-1" 


-5.4028 X 10-* 


-9.0942 X lO-i-i 


0.07 


1.2657 X 10-6 


-2.3 X lO-'^ 


-4.2719 X 10-1" 


9.3035 X 10-* 


-3.79 X lO-" 


0.075 


1.2462 X 10-6 


-2.3 X lO-'^ 


-3.9083 X 10-1" 


4.6026 X 10-* 


-4.3997 X 10-1" 


0.08 


1.2495 X 10-6 


-2.3 X lO-'^ 


-4.0554 X 10-1" 


-3.0652 X 10-* 


-3.1951 X 10-1" 


0.085 


1.2525 X 10-6 


-2.3 X lO-'^ 


-4.1011 X 10-1" 


1.0153 X lO-'^ 


-1.8389 X 10"" 


0.09 


1.2543 X 10-6 


-2.3 X lO-'^ 


-4.3337 X 10-1" 


5.0976 X lO-'^ 


-1.1825 X lO-i" 


0.095 


1.256 X 10-6 


-2.3 X lO-'^ 


-4.7088 X 10-1" 


-1.1626 X 10-6 


-6.8344 X lO-i"^ 


0.1 


1.246 X 10-6 


-2.3 X lO-'^ 


-5.1666 X 10-1" 


2 X lO--"^ 


-1 X 10-16 



Data obtained with voltages in the range of 50 — 100 mV are analyzed following the 
same procedure. The fit parameters obtained for these data are tabulated in Table [5?T1 
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As seen from the table, the coefficient 6 is a constant independent of the voltage on 
the lens as is to be expected. Its value of —2.3 x 10~^ compares well with the value of 
—2.46 X 10~^ estimated from the time period of the pendulum in the absence of the 
lens. A plot of the coefficient, c, which describes the voltage and charge dependent 
forces as a function of voltage applied to the lens is shown in Fig. I5.12[ The figure 
also contains a parabola ffited to the data. The minimum in the electrostatic force is 
determined by the balance between the charge dependent and the charge independent 
forces. The figure shows that c is minimum for an applied voltage of 75 mV. Thus, 
the background electrostatic forces acting on the pendulum at these separations are 
balanced when a voltage of 75 mV is present on the lens. 
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±6.807E-9 



Fig. 5.12: A plot of the ffiting parameter c, as a function of voltage on the lens 
calculated for the large amplitude oscillations of the pendulum. 



Fall Data: 



The other kind of data that was acquired is the data as the pendulum fall on the lens 
from separations of < 15 fim. The electrostatic forces and the Casimir force together 
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determine the motion of the pendulum in this region. 



A sample plot of this data is shown in Fig. 15. 9[ The angular acceleration is estimated 
by fitting a second order polynomial as before to every 10 points in the (t — 6) data. 
The {6 — 6) are given by Eqn. 15.441 Here too the data from multiple falls is combined 
by ordering them increasing in 9. Data is averaged over a 6 range of 1 x 10^^ radians. 
Fig. 15.131 shows a plot of {6 — 6) for the fall data with a voltage of 80 mV on the 
lens. The 'asterisks' mark the data points obtained as described above. This data is 
fitted with Eqn. I5.45[ The coefficients a, 6, c and d correspond to the various forces 
acting on the pendulum as mentioned earlier. The presence of parameters Or and d 
is essential for this fit. In the absence of d, the curvature in the data is not matched 
by the form of the fitting curve. Thus this data has signature of the Casimir force. 




Fig. 5.13: 



1.0x10 



2.0x10 

e(radians) 



3.0x10 



4.0x10 



Data: 9" Data with 80mV on the lens 
Model: 9" = a + b*9 + c/(9 + 9^) + d/9^ 
Weighting: y No weighting 



Chi'^2/DoF 



: 1.2746E-14 





= 0.99588 




a 


3.3882E-6 


±2.5474E-6 


b 


-0.00245 


±0.00307 


c 


-1.2802E-9 


±9.0135E-10 


thetaO 


0.00013 


±0.00007 


d 


-1.3015E-15 


±3.7535E-16 



{9 — 9) for fall data taken with a voltage of 80 mV on the lens and a fit to 
the data. 



The electrostatic and fibre forces have to be subtracted from this data in order to 
study Casimir force. For this, the angular accelerations obtained using the fitting 
parameters 9{fit), are subtracted from the observed accelerations, 9{obs) to get the 
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Fig. 5.14: Plot of the residual R{6)i, for the data with a voltage of 80 mV on the lens. 

residuals, i?(6') as a function of 9. Thus, for each 9i in the data set, R{9)i = 9{obs)i — 
9{fit)i. A plot of these residuals is shown in Fig. 15.141 . To minimize errors, R{9) are 
averaged over a 9 range of 5 x 10~^ radians to get average residuals and error in the 
residuals as function of average 9 as follows. 



{Oh 



2 1 2 



n 



(5.46) 
(5.47) 

(5.48) 



where n is the number of points with the 9 range of 5 x 10 ^ radians. 



The average 1/9"^ term of the fit has to be added to the {R{9))i to get the remanent 
angular accelerations with only the electrostatic and fibre background forces removed, 

9,{{9)^). Thus, 9,{{9),) = {R{9)), + {^) = {R{9)), + ^^, -, where ^™„,, and 

are the minimum and maximum value respectively of 9 in the range i that was aver- 
aged. Fig 15.151 shows a plot of remanent angular acceleration, 9i{{9)i) as a function 
of {9) with a voltage of 80 mV on the lens. ARi, represents the error bars on the data. 
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Fig. 5.15: Plot of 9i{{9)), for the data with a voltage of 80 mV on the lens. 



The estimated {6, 6) data, obtained with voltages Va in the range of 70—100 mV on the 
lens were analyzed following the same procedure to get the background subtracted, 
9i{{9),Va) and ARiiVa). The fit parameters obtained are tabulated in Table ESI 
The average value of the coefficients obtained in the fit and the errors in them are 
tabulated in Table [^751 A plot of 9i{{9i)) as a function of {9)i for three such data sets 
is shown in Fig. 15.161 A weighted average of the remanent acceleration from various 
voltages is obtained using. 



Average Remanent Acceleration, {9{{9)i)) 
A plot of this is shown in Fig. 15.171 



A Rf (alpha) 



A Fif {alpha) 

{A Rj{alphaY'^ 



(5.49) 
(5.50) 



5.4 Discussion of results 

The reduced data in Fig. 15. 171 is the force of interaction between the spherical surface 
and the fiat disc after subtracting the electrostatic forces, fibre force and constant 
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Tab. 5.2: Table of coefficients of the fit to ^ — for fall data obtained with various 
voltages applied to the lens 



Voltage on lens 
(volts) 


Fit parameters 


a 


b 


c 


0r 


d 


0.07 


2.5782 X 10-*^ 


-1.13 X 10-'' 


-1.1153 X 10'^ 


1.0 X 10-4 


-8 X 10-16 


0.075 


2.7099 X 10-*^ 


-1.59 X 10-^ 


-1.0454 X 10"^ 


8.0 X 10-5 


-8 X 10-16 


0.08 


3.3882 X 10-*^ 


-2.45 X 10-'' 


-1.2802 X 10-9 


1.3 X 10-1 


-1.3015 X 10-15 


0.085 


2.1666 X 10-^ 


-1.09 X lO-'' 


-8.2853 X 10-1" 


8.0 X 10-5 


-1.0614 X 10-15 


0.09 


2.9664 X 10-^ 


-1.9 X 10"'' 


-1.3116 X 10-*^ 


1.2 X 10-4 


-8 X 10-16 


0.095 


2.962 X 10-^ 


-2.1 X lO-'' 


-1.1186 X 10-9 


1.3 X 10-4 


-1.0701 X 10-15 


0.1 


1.3246 X 10"*^ 


-2.6 X lO-'' 


-5.7839 X 10-1° 


7.35 X 10-5 


-9 X 10-16 
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Fig. 5.16: Plot of 9i{{9)), for the data with a voltages of 75 m"V, 85 mV and 95 mV 
on the lens. 



5. Casimir Force Measurement - Strategy, Data Acquisition and Analysis 



83 



Tab. 5.3: Table of average coefficients of tlie fit to ^ — ^ for fall data obtained with 
various voltages applied to the lens 



Fit parameter 


Average value 


Error in estimation 


a 


2.585 X 10-^ 


2.541 X 10-^ 


b 


-0.00184 


2.2625 X 10"^ 




-1.011 X Kr" 


8.87796 X K)-"- 




1.01929 X 10-^ 


9.35569 X 10-6 


d 


-9.61857 X 10-16 


7.22763 X 10-1^ 




-2.5x10 



5.0x10"° 1.0x10"* 1.5x10"* 2.0x10'' 2.5x10"" 3.0x10"^ 3.5x10"* 4.0x10"* 

<6> (radians) 



Fig. 5.17: Plot {0{{9)i)) estimated from the fall data, is sho"wn. 
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Fig. 5.18: Plot of force between the lens and the pendulum disc estimated from the fall 
data(stars) is shown along with theoretically estimated Casimir force. The 
solid lines represents the Casimir force estimated assuming A = 3.26 fim 
and the dashed line represents the zero temperature Casimir force. 

(distance independent) forces and offset. This measured force plotted along with the- 
oretically expected Casimir force is shown in Fig. 15.181 The 'asterisks' represent the 
remanent force obtained after the subtraction of the electrostatic and the fibre back- 
ground forces as described in the previous section. The Casimir force is estimated 
using Eqns. 15.181 15.191 with A = 3.26 fim consistent with the theoretical estimates of 
Genet, Lambrecht and Reynaud [17] . This is averaged over a range of 5 x 10"^ radians 
in separation 6 and overploted (solid line). The Casimir force between the lens and 
the parallel plate without including finite temperature corrections is also averaged 
and plotted (dashed line). The error in the measurements at separations > 4 /im are 
less as compared to that at closer separations. This is because the pendulum spends 
less time in the region close to lens due to larger 'fall' velocity, and hence there are less 
number of data points in this region when 6 is recorded as a function of t. There is an 
offset of the measured force from the theoretically expected force even in the region 
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where the errors are smalL This offset is a small fraction of the error on the constant 
term 'a' in the fit, Eqn J5.45l to the entire data. If the data is normalized to the 
finite temperature Casimir force at 9 fim by adding a constant offset to all the data 
points, the measured force agrees with the theoretically calculated finite temperature 
Casimir force within the statistical errors of the data. The required offset of 8.4 x 10~^ 
(radians/sec)^ or 3.61 x 10"'^ dynes is ~ 1/30*'^ of the error in the estimation of the 
constant coefficient from the fit, which is 2.541 x 10~^ (radians/sec)^ [Table \5A\ . Also 
this offset is of the order of the expected residual gravitational force on the pendulum 
due to the lens [see §5.1.1]. The deviation of the data from the finite temperature 
Casimir force around separations of 4 fim is due to the fact that we have imposed a 
strict cut-off for the change over from zero temperature to finite temperature theory 
at A = 3.26 /xm. The actual force will gradually change from the dependance of 1/2;^ 
to around z = 3.26 fim. Also there will be about 10% reduction in the estimated 
Casimir force through out the range due to finite conductivity corrections. 

Thus, the measured force is consonant with the presence of a force that is ~ 1/^;^, 
viz., finite temperature Casimir force. 
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6. COMPARISON OF EXPERIMENTS ON 
CASIMIR FORCE 



Abstract: A brief review of various experiments performed to study the Casimir effect 
will be presented. We then compare the results from the recent experiments with our 
results. The first ever detection and measurement of finite temperature corrections to 
Casimir force in our experiment is highlighted. 

6.1 Experiments to Study Casimir force 

Experimental studies on Casimir force started half a century ago, a few years after 
Casimir predicted the existence of the force in 1948 [3J. A historical overview of 
the experiments was presented in the first chapter. Extensive reviews may be found 
in [U [5]. We now discuss the details about some of the pioneering experiments mea- 
suring Casimir force. 

Overbeek and Sparnaay [T5] were the first to attempt to measure the Casimir force. 
They tried to measure the force between two parallel polished fiat glass plates with 
a surface area of 1 cm^, in the distance range of 0.6 fim to 1.5 /xm. The distance 
was sometimes made 0.2 fim. The force was measured using a spring balance. The 
displacement of the spring was measured by means of a capacitive method. The mea- 
surements at 1.2 fim, 'pointed to the existence of a force which were of the expected 
order of magnitude' pLTj- The hygroscopic behaviour of the glass (quartz) surfaces 
and presence of dust particles presented difficulties. 
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Results in the distance range 0.1 /im to 1.0 /xm with several different materials were 
first obtained by Derjaguin and Abrikossova [HIT]. They measured the force between 
plano-convex fused silica lens and a fiat attached to a knife-edge balance. The deflec- 
tion of the balance was observed optically. A current in proportion to the deflection 
of the balance from its mean position was generated and passed though a coil. The 
coil was pivoted between the pole pieces of a magnet. The magnet and the coil were 
mounted such that the current flowing through the coil produced a turning force on 
the balance opposite to that due to attractive molecular forces and kept the balance 
from turning. This current was, thus, a measure of the force. The separation be- 
tween the lens and a fiat was measured using "Newton's rings". They could show 
that Casimir's formula provided a better fit to the experimental data in contrast 
to London's theory for the van der Waals forces. They also did experiments with 
two quartz plates, quartz and chromium plates and with crystals made of thallium- 
bromide and thallium iodide. It was these experiments that, in a way, motivated 
Lifshitz to derive a macroscopic theory of molecular forces [1]. 

Sparnaay [16] repeated his measurements with metal plates in 1957. He measured 
the force between chromium plates and chromium steel plates. The measured scheme 
was similar to the previous experiment with a spring balance of sensitivity between 
(0.1 — 1) X 10^"^ dynes. The distance was varied between 0.3 /im and 2 fim. The 
measurements did not 'contradict' the expected force per unit area from Casimir's 
relation. But large systematic errors and electrostatic forces prevented a detailed 
quantitative study. 

The van der Waals force between two curved mica surfaces were measured in the dis- 
tance range of 1.5 nm to 130 nm by Israelachvili and Tabor [101 E]. They measured 
the force between two crossed cylindrical sheets of mica with the effective contact 
resembling that between a sphere and a plate. The surfaces were silvered and the 
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distance between them was measured by observing interference patterns with white 
hght. The force was measured by two methods. In the short distance range of 1.5 nm 
-20 nm a 'jump method' was used. One surface was rigidly fixed to a moveable base 
and the other was mounted on a cantilever spring. When the fixed surface was moved 
towards the cantilever, at some point - depending on the stiffness of the spring - the 
two surfaces jumped into contact. The measurement of this jump distance as a func- 
tion of spring distance was used to measure the force of attraction. For the distance 
range of 10 nm - 130 nm a resonant method was used. One surface was supported 
on a rigid piezo-electric crystal and was set vibrating at very small amplitudes over a 
convenient range of frequencies. The other was supported facing it on a stiff spring. 
The natural frequency of the latter depended both on the spring constant and on the 
van der Waals forces exerted on it by the other surface. By determining the resonant 
frequency as a function of separation, the force law was deduced. They showed that 
there is a gradual transition between the non-retarded and retarded van der Waals 
as the separation is increased from 12 nm to 50 nm. 

The next major set of improved measurements with metallic surfaces were performed 
by van Blokland and Overbeek [18J in 1978. They measured the forces between a lens 
and a flat plate coated with chromium using a spring balance at distances between 
0.13 fim and 0.67 /im. Precautions were taken to characterize and eliminate electro- 
static forces by careful experimentation. Precise measurement of the separation was 
made by measuring the lens-plate capacitance. A detailed comparison of the data 
with Lifshitz theory, taking into account finite conductivity effects was done before 
concluding that "the measured force and the calculated force are in excellent agree- 
ment". It can be considered as the first unambiguous demonstration of the Casimir 
force between metallic surfaces. 

After couple of dormant decades, interest in the measurement of forces in the sub- 
millimeter range was aroused by theories that predicted the existence of new forces 
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in this range. The earhest of these experiments was by Lamoreaux in 1997 [21] and 
around the same time we started the design of our experiment [S]. The apparatus 
built by Lamoreaux consists of a glass plate coated with gold over copper suspended 
as a torsional pendulum. Its attraction towards a spherical lens also coated with 
copper and gold was measured for distances between 0.6 /xm - 6 /im. The plate was 
held fixed with respect to the spherical surface by a voltage applied to two compen- 
sating capacitor plates held parallel to the 'Casimir' plate. By measuring the change 
in the voltage required to hold the plates parallel when the lens was moved towards 
the 'Casimir' plate, the Casimir force was measured. He detected and corrected for a 
contact potential between the lens and the plate of 250 mV. The measured Casimir 
force agreed with theory at the level of 5% up to separation of about 2 fim. The sen- 
sitivity of the experiment was not adequate to measure the force beyond this range. 

In a later experiment, Mohideen and Roy [H] measured the Casimir force between 
a metallized sphere of diameter 196 /im and a flat plate for separations from 0.1 to 
0.9 fim. They used the high sensitivity of an Atomic force microscope to report a 
statistical precision of 1% at the smallest separations. The electrostatic force due a 
residual voltage of 29 mV was measured and subtracted from the data. 

A micromechanical actuator based on Casimir force was developed by researchers at 
Bell Labs [4] in 2001. The device, fabricated using standard nanofabrication tech- 
niques, consists of a 3.5 fim thick, 500 /im square heavily doped polysilicon plate 
freely suspended on two opposite sides by thin torsional rods. The other ends of the 
torsional rods were anchored to the substrate. There was a 2 fim gap between the 
top plate and the underlying substrate. A gold coated ball was suspended above one 
side of the plate, which caused the plate to tilt due to Casimir force when the plate 
was brought close to the ball. This tilt was determined by measuring the capacitance 
between the plate and the substrate accurately. A residual voltage of 30 mV was 
determined and compensated for. The measurements matched with Casimir force 
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theory at the 1% level. The errors were mainly due to the corrections for finite con- 
ductivity and surface roughness that have to be applied to the theory. 

Recently, Bressi et al. p] have measured the Casimir force between parallel metallic 
surfaces. This is the first and only precision measurement of the Casimir force in 
the 'parallel plate' configuration. The force was exerted between a silicon cantilever 
coated with chromium and a similar rigid surface and was detected by looking at 
the shifts induced in the cantilever frequency when the latter was approached. The 
motion of the cantilever was monitored by means of a fiber optic interferometer. Elec- 
trostatic force due a residual voltage of ~ —68 mV is systematically determined and 
corrected for. The scaling of the Casimir force with the distance between the surfaces 
was tested in the 0.5 /im - 3.0 /im range, and the related force coefficient was deter- 
mined at the 15% precision level. 

Very recently, the Casimir force between two dissimilar metals was measured by Decca 
et al. |6j in the separation range of 0.52 /xm - 2.0 yum. They measure the force of 
attraction between a copper layer evaporated on a micro-mechanical torsional pendu- 
lum and a gold layer deposited on an Aluminium oxide sphere with 600 /im nominal 
diameter. The force is inferred from the deflection angle of the torsional oscillator 
determined by measuring the capacitance between the oscillator and two fixed elec- 
trodes. The sensitivity of the force measurement was further improved by measuring 
the change in the resonant frequency of the oscillator due to the presence of the 
sphere. Residual voltage difference of 632.5 mV was determined by experimentation 
and corresponds to the difference in work functions of the gold and copper layers. 
The experiment highlights the need for simultaneous measurement of the dielectric 
constant of the metallic films used for a better understanding of the Casimir force 
between non-ideal bodies at these separations. 
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6.2 Comparison of results 



The Casimir force has been measured over a wide range of separation from 1.5 nm to 
3 /im by various techniques. Most experiments have been performed in the sphere- 
plate geometry rather than in the plate-plate geometry due to stringent conditions of 
parallelism required in the latter. All the experiments encountered a background of 
strong electrostatic forces due to residual voltages present on the interacting bodies 
even when these are well grounded. The Casimir force for the sphere-plate geometry, 
depends on the radius of the sphere used and is typically much larger than the electro- 
static force at sub-micron separations. But at larger separations, the measurements 
are dominated by the electrostatic forces, making it difficult to perform experiments 
in the several microns range. For example, an electrostatic contribution comparable 
to the Casimir force at 1 fim becomes 10 times or more stronger than the finite tem- 
perature Casimir force at 10 fim. The experiments performed thus far, have been 
able to detect effects due to finite conductivity and surface roughness of the metals 
on the Casimir force but the finite temperature correction has not been observed. 

Our experiment measures the Casimir force in the separation range of 2 fim to 9 fim. 
This region extends well into spacings beyond ~ 3.3 /im where corrections to the 
Casimir force due to the ~ 300° K environment become significant. As expected, 
the data clearly indicate the presence of the effects of finite temperature. Data from 
recent experiments that measured the force in the sphere-plate configuration are sum- 
marized in Fig. 16. 1[ The figure shows the force measured in these experiments divided 
by 2'kR, where R is the radius of curvature of sphere used in the respective experi- 
ments. The results presented in this thesis are represented by the 'stars'. The data 
from various experiments overlap in the region where the errors are small. Our data, 
below 3.3 fim overlaps well with the measurements by Lamoreaux. The slope of the 
data represents the power of the force law obeyed by the data. Around separation 
of ~ 4 fim the plot distinctly shows the change in slope from ~ —3 to ~ —2 within 
the errors of the experiments. There is a clear indication of the change-over in the 
Casimir force law from the zero temperature theory to that due to finite temperature 
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Fig. 6.1: Data from recent Casimir force experiments using the sphere-plate geometry 
by Lamoreaux [21] , Mohideen and Roy [T^ , Chan et. al. [1] , Decca et. al. [6] 
and data from our experiment. 
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around this separation. 

In summary, our current study extends the observation of the Casimir effect up to 
separations of 9 /im, i.e., well into the region where finite temperature effects become 
important. The experiment confirms the existence of these effects and agrees with 
the standard theory at the 20% level. 

Keeping in mind the wide-ranging importance of Casimir forces it would be useful to 
conduct such studies with greater accuracy and also extend them to longer separations 
of ~ 50 //m or more. 
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7. BOUNDS ON THE STRENGTH OF NEW 
MACROSCOPIC FORCES AND FUTURE 
DIRECTIONS 



Abstract: We briefly review the experimental bounds on the strength of inverse square 
law violating short-range interactions. The contribution of our experiment in this 
direction will be presented. The future scope of our experiment to test and constrain 
the predictions of fundamental theories will be discussed. 



7.1 Constraints on new macroscopic forces 

Extensions of the Standard Model predict the existence of a variety of neutral light 
bosons. The exchange of these will mediate forces that lead to deviations from the 
inverse square law of gravity at distances related to the mass of the bosons. String 
and M-theories that attempt to unify the fundamental forces close to the Planck scale 
or the theories with large extra dimensions that attempt to overcome the "hierarchy 
problem" by invoking unification of the fundamental forces at the TeV scale of electro- 
weak symmetry breaking, predict variations from the inverse square law of gravity at 
sub-millimeter distances. An overview of these effects was presented in Chapter 1. In 
general, the new macroscopic forces can be described by an addition of a Yukawa type 
potential term to the gravitational interaction between two point masses as given by, 

= + (7.1) 
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where a represents the couphng strength of the interaction and A the range. A recent 
review of the experimental and the theoretical status of the inverse square law tests 
can be found in [jlj. Given this form for the potential, constraints can be placed on 
the parameter space of a — A from experiments that study long-range interactions 

7.1.1 Astrophysical Bounds 

Astronomical tests provide the best constraint on the parameter a at distance scales 
A > 1 km. These are typically obtained from observation of the Keplerian orbits of 
the planets and the satellites and looking for deviations from normal gravity. The 
measurements of precession of the lunar orbit using Lunar Laser Ranging studies 
provides very good constraints in the 10^ m scale. Precession is expected due the 
quadrupole field of the Earth, gravitational perturbations from the other planets in 
the solar system and general relativistic effects. Presence of Yukawa type interactions, 
will add to the precession and hence limits can be placed on a after accounting for the 
expected sources of the precession. A summary of the constraints from astrophysical 
data is reproduced from [1], [15] in Fig. 17.11 In the range A < 1 m, these limits have 
been surpassed by laboratory experiments [22| [25] . Details on geophysical and 
laboratory constraints in the range A > 1 m can be found in [77] and [5]-|llj. 

Astronomical data can also provide evidence for the presence of extra dimensions 
[2]. The cooling rate of the supernovae will increase by the presence of the extra 
dimensions because energy from the explosion will be radiated into these dimensions 
as well. This would lead to the reduction in the number of neutrinos emitted from the 
supernovae. SN 1987A data constrains the size of the extra dimensions to < 0.7 /xm 
when the number of extra dimensions is 2 [2l [121 ISOl [21]. This does not imply 
that experiments at sub-millimeter scales will not observe effects due to the new 
dimensions. A single large dimension of size 1 mm with several much smaller extra 
dimensions is still allowed. 
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Fig. 7.1: Constraints on inverse square law- violating Yukawa interactions with A > 
1cm. The Lunar Laser Ranging (LLR) constraint is based on the anomalous 
precession lunar orbit; the remaining constraints are based on Keplerian 
tests. This plot is reproduced from pTj. The region in the {a, A) plane above 
each curve is excluded, and below each curve is allowed. 

7.1.2 Bounds from Laboratory tests of Inverse square law 

Constraints on a in the sub-millimeter to centimeter scales can be obtained from lab- 
oratory experiments of the Eotvos and Cavendish type. The Eotvos-type experiments 
test the equivalence principle by measuring the acceleration imparted by Sun, Earth 
or some laboratory attractor to various materials of the same mass. The presence 
of any additional force that couples to the material properties other than the mass 
will show up as a difference in the acceleration experienced by the various materials. 
The Cavendish-type experiments are direct tests of the Newton's inverse square law. 
The tightest constraints in the 1 mm to 200 fim region to date are provided by such 
experiments. These are shown in Fig. 17.21 (reproduced from [T]). The shaded region 
in the (a. A) plane above the continuous curves is excluded. 
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Fig. 7.2: Constraints on inverse square law- violating Yukawa interactions with 1 /xm 



< A < 1 cm. This plot is from [T]. The heavy curves give experimental 
limits. The curve marked Irvine is a Eotvos type experiment. 'Lamoreaux' 
is a Casimir force measurement. The other curves are from Canvendish-type 
experiments. The theoretical expectations from various scenarios are also 
shown. 



7.1.3 Bounds from Casimir Force Measurements 

For distances of < 0.1 mm, Casimir force measurements provide the tightest con- 
straints. Casimir force has generated extensive experimental and theoretical interest 
in the last decade. Precision experiments have been performed in the distance range 
of about 0.1 fim to 3 fim and the theoretical correction to the Casimir force due to 
finite temperature, finite conductivity and surface roughness have been calculated 
accurately. This has enabled comparison between theory and experiment at the 1% 
level for the smallest separations. If the inverse square law violating interactions are 
present due to any of the scenarios mentioned in Chapter 1, they will show up as 
additional forces in the measurement. Thus, constraints on the strength of these 
interactions can be obtained by looking at the deviation of the measured Casimir 
force from the theoretically expected value and attributing it to Yukawa type inter- 
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actions [31 HJ [lU [I3]. A plot of the constraints to date is summarized in Fig. I7.3[ 
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Fig. 7.3: Plot of the constraints on the Yukawa interaction parameter a from various 
ranges of A reproduced from Curves 7-10, 12 follow from Casimir force 
measurements, Curve 11 from van der Waals force measurements. Curve 
6 is from an experiment that measured deviations from Newton's law of 
gravity [2S]- The typical prediction of extra dimensional physics is shown in 
Curve 13. The region, in the (a. A) plane, above each curve is excluded and 
the region below each curve is allowed. 

7.1.4 Estimation of bound from our experiment 

Our experiment measures the Casimir force between a spherical lens and a flat disc 
made of 'BK7' glass (density, = 2.51 g.cm~^), coated with gold layer of thickness, 
6 = 1 yum (with density, ps = 19.32 g.cm"'^). The radius of curvature of the lens, 
-R = 38 cm and its total height Xmax = 0.07 cm. The radius of the flat disc, / = 4 cm 
and its thickness, g = 0.4 cm. The force due to Yukawa type potential for this 
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configuration is given by (see Appendix C for details), 



Ta{z) = ATc^GaRX^ [p., - {ps - Pv)e~'/^] Ps 6 e-'/\ for A « 6, (7.2) 
J-fe(^) = Arc^GaRX^ [p., - {ps - P.)e"'/^] ( p. 5 + p. A e~'/^ ) e~'/\ (7.3) 



for A >> 6. 

where a and A are parameters of the theory as defined by Eqn. 17.11 

The upper limit of constraints on the Yukawa-type interactions from our experiments 
can be estimated by two methods. A conservative estimate of the parameters of the 
Yukawa-type interactions can be obtained by assuming that all of the force, fobs{z) 
measured in our experiment is due to the Yukawa-type potential. Thus, 



We solve Eqn. 17.41 for the observed values of the force at z = 1.82 pm, 3.08 pm, 
4.22 pm, 5.36 pm, 6.51 pm, 7.641 pm, 8.77 pm to get a as a function of A, we call 
this aconservative{X) . Siuce the thickness of the gold coatings, 6 = 1 pm, we derive 
ctconservativeiX) in the range A = 2 pm to 16 pm using Eqn. 17.31 This gives us seven 
possible functions for aconservative{X) corresponding to the seven separations at which 
we have measurements. A plot of these is shown in Fig. 17.41 The lowest values of 

(^conservative 

(A) are obtained from data at 2; = 8.77 pm in the range 7 pm < A < 16 pm 
and from data aX z = 4.22 pm in the range 3 pm < X < 7 pm. 

An optimistic estimate of the parameters of the Yukawa-type interactions can be 
obtained from the error in the measurements of the force. Our experimental data 
matches with the finite temperature Casimir force between the lens and the plate 
within the errors of force measurement at separation z > 4 pm. The dominant con- 
tribution to the finite temperature Casimir force for our configurations [see Appendix 
A], at these separations is given by. 



fobs{z) 



(7.4) 




for z > Xt = 3.26 pm; 



(7.5) 
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Fig. 7.4: aconservativeW estimated from the total force measured in the experiment. 

The corrections to Casimir force due to surface roughness are < 0.6% in the 4.5 fim to 
9 /im range [3] and can be ignored. The finite conductivity of gold surface gives rise 
to a reduction in the Casimir force by 10% over the entire range of the experiment 
[T71 [T6] . Thus, the net theoretical force between the lens and the plate is given by, 

fth = fciiz) + ApUz) + J^{z); (7.7) 

where Apfci{z) represents the corrections due to finite conductivity and JF is the force 
due to hypothetical Yukawa-type interactions. This is derived in Appendix C for our 
geometry and is given by Eqns. 17.21 and 17.31 Thus, the Yukawa parameters can be 
estimated using, 

\f,^(z)-Uz)\ ^AFiz); (7.8) 

where AF{z) is the error in the measurement. If absolute error in our force measure- 
ment at a distance z is e{z), we can solve Eqn. 17.81 assuming AF{z) = 3e{z). The 
corrections due to surface roughness become less important for larger separations and 
for z < 4 fim, the theoretical Casimir force is not strictly defined by Eqn. 17.51 We use 




7. Bounds on the Strength of New Macroscopic Forces and Future Directions 



103 



z = 4.22|j 
z = 5.36|j 
z = 6.51n 
z = 7.64n 
z = 8.77n 



-5.4 -5.2 -5.0 -4.8 

Log (W1m) 

Fig. 7.5: aoptimisticl-^) estimated from the errors in the measurement. 

the data at 2; = 4.22 /im, 5.36 fim, 6.51 fim, 7.64 fim, 8.77 /im and solve for a from 
Eqns. 17.71 and 17.81 in the range 3 /im < A < 16 fim. The values of a obtained by this 
method are represented by aopUmisUc and are plotted in Fig. 17.51 The lowest values of 

(^optimistic 

(A) are obtained from data at 2 = 8.77 /im in the range 6 /zm < A < 16 /xm 
and from data ai z = 7.64 /^m in the range 3 /im < A < 6 /im. 

In Fig. El, the lowest values of aconservativeW and aopHmisticW estimated as de- 
scribed above are plotted along with those from other experiments as presented by 
Mostepanenko in 

The lines marked 'Zerolab' are constraints as calculated from data presented in this 
thesis. Curve 6 is from the Canvendish-type experiment [25]. Curve 7, a is from the 
torsion balance experiment to measure Casimir force [21]. This is an upper limit of the 
constraint we can obtain from our experiment. The constraints from the conservative 
estimate are already at the level of those obtained from [24J. The optimistic estimates 
give the lowest values of a obtained so far in the range 3 /xm < A < 6 /xm and also the 
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Fig. 7.6: Constraints on the inverse square law violating Yukawa interactions obtained 
from data presented in the thesis ('Zerolab') along with the best constraints 
in the same range of A from [26] 

best ever constraints from Casimir force measurements in the range 3 /im to 16 /im. 
By calculating the theoretical forces precisely, the constraints can be improved and 
also extended to shorter distances. This also means that by choosing the materials 
of the pendulum and attracting plate appropriately (for example, thicker coatings, or 
higher density materials) we may be able to improve the present constraints by more 
than an order of magnitude in the distance range 3 fim to 10 /im. 
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7.2 Future Directions 

We are in the process of performing a next generation experiment, with the same gen- 
eral configuration but improved precision by reducing the systematic and statistical 
noise. In the new scheme we will perform a null experiment where the deflection of 
the pendulum due to the torque from the lens is balanced by the capacitive torques 
from a pair of metal plates located close to the disc of the pendulum (see §2.3). The 
voltages on capacitor plates are controlled through a servo-loop and the servo-signal 
is a direct measure of the torques acting on the balance. This mode of experiment 
has several advantages. The gain in the feed-back helps to reduce the noise in the 
system (see for example [ISj)- The statistical accuracy can be considerably improved 
by increasing the integration time. The errors associated with the estimation of ac- 
celeration from the limited time series data are avoided. The experiment can be 
performed at closer separations as well. The torque due the lens can be modulated 
by moving the lens back and forth, further reducing the noise in the bandwidth of the 
observation. Efforts are also on to characterize the background electrostatic forces 
better and hence decrease systematic errors. 

With the lens-plate combination, Casimir force at separations greater than 30 fim 
will be difficult to measure as electrostatic and gravitational forces will start to dom- 
inate. A similar experiment would be repeated with a fiat plate instead of a lens for 
separation up to about 100 /im. The parallelism between this plate and suspension 
disc is less critical at these separations. The dependance of the Casimir force on 
separation will be as opposed to the typical l/d"^ or 1/d dependance of the other 
background forces. Such an experiment will be able to put stronger constraints on 
the new forces with strength close to that of gravity in the sub-millimeter range as 
the parallel plate configuration maximizes the sensitivity to these forces as opposed 
to the plate-sphere experiments. 

We expect a 5% accurate measurement of the Casimir force in the range 10 /xm to 
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100 //m, and also improved constraints on hypothetical Yukawa-type forces with range 
10 /im to 100 /im, from these future experiments. 
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APPENDIX 



A. CASIMIR FORCE BETWEEN INFINITE 
PARALLEL PLATES 



The Casimir force is a manifestation of zero point energy of the electromagnetic field. 
In the quantum mechanical description of the electromagnetic field, the allowed energy 
levels of the electromagnetic wave of angular frequency uj are given by the Planck 
relation En = {n + ^)hu, n = 0, 1, 2, 3 ... . The integer n, for the electromagnetic 
field, corresponds to the number of photons. The n = or the 'vacuum' state also 
has an energy of l/2huj associated with it. Thus, vacuum is not empty and contains 
fluctuations of the electromagnetic field. 

Let us now calculate the effect of two parallel conducting plates separated by a dis- 
tance d, placed in such a field using the 'mode counting' technique [Sj. The field 
modes within this cavity (Fig. lA.ip are different from the free space modes. In ad- 
dition to satisfying the wave equation, V^A = ^^r, together with, V.A = 0, 
the field inside should have the tangential component of electric field vanishing at the 



z= z = d 

Fig. A.l: Two perfectly conducting infinite plates, placed distance d apart. 
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perfectly conducting boundaries. The solution that satisfies these conditions is given 

by, 

A{r) = A,{r)x + Ay{r)y + A,{r)z; (A.l) 
/ g 

where Ax{r) = y y o-x cosk^x sinkyy sink^z, (A. 2) 

/ g 

^3/(^) ~ y V '^^ sinkxX cos kyy sin/c^z, (A. 3) 
/ g 

Az{r) = y Y o-x sink^x sinkyy coskzZ; (A. 4) 

with, 

al + al + al = l;V = L'^L, & (A.5) 

111 7im 7m 

kx = — ; ky = k^ = — ; where l,m,n = 0, 1, 2, 3, . . . (A. 6) 

V.A = also requires that 

k^Ax + kyAy + k,A, = or (A.7) 
^{lA, + mAx) + ^{nA,) = 0. (A.8) 

There are two polarizations possible for each value of l,m,n unless one of them is 
zero. In that case, only one polarization is possible due to Eqn. IA.8I Thus, the 
allowed modes within the cavity are. 



^Imn = hmnC = ( — + — + — ) . (A.9) 



2 



The zero-point energy of the field inside is. 



E{d) = Y; 2^nwimn = E' + ^ + z^) ' • (^-^o) 

Imn Imn ^ ^ 



The factor 2 arises from the 2 independent polarizations for l,m,n ^ and the prime 
on ^ indicates that when / or m or n = 0, a factor of half has to be inserted. 
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When L >> = d. 
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(A.ll) 
(A.12) 



This is infinite. Outside the cavity, will also be replaced hy ^ J dkz and 



TT J TT , 



OO poo /*oo 



(A.13) 



'0 ^0 ^0 

This is also infinite. But the potential energy of the system when the plates are 
separated by a distance d is U{d) = E{d) — E{oo). 

2 



U{d) 



L 



TT 



he 



E' 



.n=0 

^ poo poo poo 
TT 



^0 



2 2 n'^ir'^^ ^ 
dkxdky ^ k^ -\- ky -\ 



dk^dkydk^ ij^x 



(A.14) 



'0 ^0 ^0 

In polar co-ordinates, m, ^ in kx^ky plane, dkxdky ~ ududO, 9 ranges from to 7r/2 
for positive kx, ky, and we have. 



U{d) 



nj 2 

d 



TT 







E'X 



(A.15) 



This is still infinite. Wc now introduce a cut-off function f(k) = f{[u^ + /c^]^/^) 
such that f{k) — 1 for k « km and f{k) — for k » km- We could say that for 
wavelengths comparable to atomic dimensions, the assumption of perfect conductivity 
breaks down and hence a cut-off is necessary {k^. ~ I/oq, is Bohr radius). Thus, 
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Eqn. IA.15I becomes, 
U{d) 



d 

Jo 



{n=0) 



du u { + 



rf2 



dK / duu{u'^ + k'^^~^ f{[u^ + k'^ 



Defining x = u d j'n & k = hzdj-n, 



4:71 d^ 



1 TT 



Now, 



t/(ci) 

where 



-I °° /-oo 

-F(0) + ^F(n)- / 

n=l 

/ dx{x + K^y f{^[x + K^]^). 

Jo d 



r"00 



According to Euler-Maclarin summation formula [T], 



En 

n=l 



'dKF[K) = -1f(0) - ^F'(O) + ^F"{0) + 



for F(oo) ^ 0. To evaluate -F"(0), we note that. 



(A.16) 
(A. 17) 



POO POO _ 

dn dx{x + f{-[x + k'^]^^) . (A.18 
Jo Jo "J 



(A.19) 
(A.20) 



(A.21) 



poo 

FiK) = / du v^/(-v^), F'iK) = -2K'f{-K). (A.22) 
Jk2 d d 

Then F'(0) = 0, F'"(0) = -4 and ^''(O) = forn > 3, if the cut-off function vanishes 

at K = 0. Thus, 



-(^) - ^^-^ (-4) 



720c/3 

which is finite and independent of the cut-off function. The attractive force per unit 
area, 

- (A.26) 
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^2 



z = 



z = d 



Fig. A. 2: Two semi-infinite dielectric slabs of dielectric constants ei and 62, placed 
distance d apart in a medium of dielectric constants €3. 

which is the Casimir force. Thus, Casimir showed that changes in zero-point energy 
can be finite and observable. 

A.l Force between Dielectrics 

In the previous section, we went through Casimir's derivation of the force between 
two perfectly conducting parallel plates. In experimental situations, the simplified 
assumption of perfect conductivity at all field frequencies is unrealistic. The dielectric 
properties of the media should also be included. Liftshitz developed the first 
macroscopic theory of forces between dielectrics. His results reduce to the Casimir 
force derived above for the case of perfect conductors. We will now derive Lifshitz's 
results in a manner similar to Casimir's approach of mode counting |5|. 

Consider the case of a medium with dielectric constant €3(00) sandwiched between 
two semi-infinite media of dielectric constants ei{uj) and €2(0;) [Fig. IA.2j . We will 
calculate the force per unit area from the total zero-point energy of the modes Un 
that are possible in this situation. The modes that are possible are the solutions of 
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the Maxwell's equations, 

V.^ = 0, 

^ ^ ^dB 
c at 

V.B = 0, 

c at 

where, we have assumed the media to be isotropic, with magnetic permeability fJ> — 1, 
and zero net charge density, so, D{f, t) — e{uj)E{f, t). These modes should also satisfy 
the appropriate boundary conditions. Consider solutions to the Maxwell's equations 
of the form E{f,t) = i5o(r, t)e~*'^*, B{f,t) = 5o(^, such that in each region, 

V.Eq, V.Bq = and, 

V^Eo + '4<'^)Eo = 0, (A.27) 
V'So + ^eHA) = 0. (A.28) 

Across the boundaries, 

(i) normal component of D{= eE) should be continuous, 

(ii) tangential component of E should be continuous, 

(ill) normal component oi B{— H) should be continuous, and 
(iv) tangential component of B should be continuous. 

Let us assume solutions of the form, 

= [e,(z)x + e,(^)y + e,(z)i]e'(*^-^+*^^^), (A.29) 

Bo(f) = %{z)x + hy{z)y + h,{z)z\e'^^-''^^^y\ (A.30) 

sothat,^-i^2e. = 0, ^-K%,^Q, (A.31) 

dz'^ dz^ 

where = kl + k^- —e{uj) (A.32) 
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and likewise for y and z components. In free space, e{uj) = 1 everywhere for all uj and 
the usual plane wave solutions are obtained for K"^ < 0; kl + k^ + kl = K"^ = — fc^. 

For the case under consideration, let us assume e{uj) to be real and > in all three 
media. We can chose a coordinate system in which ky = and = k"^ ~ t{uo)^ , k = 
kx- 

V-Eq = 0, implies that ike^ + -j^ = 0. (A. 33) 

V X 4) = i^Bo, gives 
i?o(r) = 

V.-Bo is also identically satisfied by Eqn. IA.34I 
The boundary conditions mentioned above imply, 

(i) e{uj)ez{z) should be continuous, 

(ii) Cx should be continuous which is satisfied if ^ is continuous (from Eqn. IA.33I) . 

and 

(iii) Cy should be continuous, 

(iv) ^ and ke^ + should be continuous. 

But, from Eqn. IA.32I and IA.331 kcz + = \^{l^)^^z- Thus, all conditions are 
satisfies if, 

a. t{ijj)ez{z) and ^ are continuous, and 

b. Cy and ^ are continuous. 

Now, since — K^Cz = 0, in each region, ignoring exponential growing solutions, 
we have, 

ez{z) = Ae^'\ for2<0; (A.35) 

= Be''''' + Ce-'''''' for < 2 < 0; (A.36) 
= De-''^'^ forz>0; (A.37) 



. c dcy ^ c 
I -^x 

CO dz CO 



kcz + ) y + —kCyZ 
dz J ijj 



Akx 



(A.34) 
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where Kj = fc^ — ej{u)^. Applying the boundary conditions 'a', we have, 



- eiA + esB + esC = 0; 

-KiA + KsB-K^C = 0; 

ege^^^E + ese-^^^C - eae^^'^Z^ = 0; 

Kse^--"^B - Kse-'^'^'^C + K2e^^'^D = 0. 



(A.38) 
(A.39) 
(A.40) 
(A.41) 



For non-trivial solutions to exist, 



This yields. 



-ei 




63 

























-626" 


-K2d 







-7^36-^3'^ 


K2e- 


K2d 




+ eiKs 


€3^2 + e2-ft'3 
€3/^2 - ei^3 


^2Kzd 


- 1 



0. 



Similarly, it can be shown that the boundary conditions 'b' are satisfied if, 

-1 = 0. 



+ K2 + K^ 2Kzd 

K, - K:,- K2 - Ks'"" 



(A.42) 



(A.43) 



Equations IA.42I and IA.43] are the conditions on the allowed modes uJn- Both cannot 
in general be satisfied simultaneously. For boundary conditions a and b to be satis- 
fied simultaneously; (a) Eqn. IA.42I should be satisfied with e^^ = or (b) Eqn. IA.43I 
should be satisfied with = 0. Thus we have two kinds of modes. Since these are 
exponentially decaying functions of z for z < and z > d, these are called surface 
modes [2]. 



Let us now calculate the zero-point energy of these modes. The energy is given by. 
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where ujna are the modes of type (a) and uJnb are the modes of type (b). Since and 
ky are continuous, the summations over them in uJnai a = a,b can be replaced with 
integrals. 

E ^ (^) V /^^^E = 1 2-kdkJ2, (A.45) 

where L is the length along x and y sides of the 'quantization box' and denotes 
sum over the solutions of [02] and [03l Thus, 

fc r2 roo 

E{d) = ^ dk. 

In order to evaluate this summation, we employ the 'argument theorem' from the 
theory of functions of complex variables: for a function f{z) which is analytic except 



E^^«(^) + ^^^Nb{k) 



(A.46) 



for poles on and inside a simple closed circle C, 



(A.47) 



where is the number of zeros and P the number of poles f{z) inside C (see for 
example, [3]). It can be shown from this that. 



27Ti 



c 



-dz 



f(Zi) = 0O 



= sum of zero's of f{zi) inside C 
— sum of poles of f{zi) inside C. 

Let Fa{u) and Fb{u) be the L.H.S of Eqns. |X42] and |X43] then, 

'y^^UNajK) = sum of zero's of Fq,, a = a, b. 

N 



(A.48) 

(A.49) 
(A.50) 



(A.51) 



The poles of Fa, considered as function of complex variable u, are independent of the 
boundaries and hence of d. Therefore, 

^— / uj ^^^\ duj = ^ ujNa{K) — (Term independant of d). (A. 52) 

The ujNa{,K) of interest are those that are positive and real. Hence, the closed curve 
C is defined by the imaginary axis of complex uj plane and a semicircle in the right 
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half of this plane, with radius extending to infinity. As the d dependant term does 
not contribute to the force, for this analysis, we can write. 



E{d) 



hL^ 1 

47r 2TTi 



dk k 



F'iu) 



ijj- 



c 



Faiuj) 



FLiu) 



c 



FJu) 



(A.53) 



Each of these contour integrals, can be written as a part along the imaginary axis 
and a part along the infinite semicircle. The later integral is c?- independent and does 
not contribute to the force. The former integral is. 



■,^_dF}M)_ 



—I 



—I 



d 



di ^-logG«(e); 
di logGaiO': 



(A.54) 
(A.55) 
(A.56) 



where Ga{C,) = Fa{i^), a = a,b. Writing exphcitly 
GbiO - 



e^Ki + eiK^ €3X2 + £2/^3 ^2K3d _ ^. 
esKi - eiKs ' €3X2 - eiK^ ' 

K1+K3 K2 + K3 2K,d _ ^. 

K^-K3K2-K3 



(A.58) 



where, now = and K] = e + ejii^^. The zero-point energy associated 



--jy-^j -J 

with the allowed modes is then, 



hL 



2 r°° 



dk k 



d^ logGaiO 



d^ logGaiO 



+ d independant term. 

(A.59) 



The force F{d) = -£E{d), 



Fid) 



hL^ d 
"S^dd 



dk k 



/oo 
di logGtiO 
-00 



(A.60) 



But, from equation |X37] and |X38l ^^G^iO = '^^3 (G^ + 1), and 



d_ 
dd 



/oo poo 
dUogG^iO = J d^—logG^iO; 

POO 1 



d^ 



00 
00 



GaiO 



[2K3 (G„ + 1)] 



d^K3 



d^ 



K3 



Ga{0' 



(A.61) 
(A.62) 
(A.63) 
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The first term is independent of the presence of the dielectric slabs and hence is not 
related to the force between them. Therefore, the force per unit area, 

n 



87r2 



/"oo rc 
/ dkk 
Jo J-c 



1 1 

+ 



Kj and are even functions of ^, and so is Ga', Hence 



F{d) 



h 



More explicitly. 



F{d) = 



h 



POO PC 

/ dkk 

Jo Jo 



diK^ 



1 1 

+ 



27r2 



dk k 



d^Ks 



-1 



+ 



K^-K3K2-K3' 



where, Sj — ej(i^) and — k'^ + ej{i^)^. 

To compare with Lifshitz theory we introduce a variable p such that. 





^2 

s / 2 


-!)■ 


Then, Kl = k^ + 


e 2 




andir?2 = + ^1,2^ 
' c 


s / 2 






e ( 2 


-1 + ^; 

^3 / 




_ e 2 

= ^3-^Sl,2, 


where 2 = 



£1^ 
£3 ' 



^2 

Changing to these set of variables, dk k = e^^dp p and. 



F{d) 



27r2c3 



dp p^ 



33/2 



,Sl -P'S2 -p' 



(A.64) 



(A.65) 



(A.66) 



(A.67) 
(A.68) 



(A.69) 



(A.70) 
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For the case when there is vacuum between the slabs, (63 = 1) 

h 



27r2c3 



^1 + ^2 + e2P ^2^pd/c _ ^ 

Si - eip' S2 - eip' 



\si-p S2-p J 

which matches with Lifshitz's result [1]. 



-1 



(A.71) 



A. 2 Casimir force at finite temperature 

When we derived the Casimir pressure between the dielectric slabs, we assumed the 
vacuum state for the field and added only the zero-point energy of the allowed modes 
to calculate the energy of the field. At any finite temperature, 

En = (n{uj) + - j fkj, where n^uj) = — . (A. 72) 

Thus to calculate the energy of the allowed modes at a finite temperature, we have 
to use Eqn. IA.72t rather than just |/k^, 

2 

Eid).2.{^ + ^1; (A.74) 



2 



= S(<i)cothJ^. (A.75) 

The easiest way to now derive the force expression is to use u = iC, in Eqn. IA.71I and 
multiply by the factor coth 2^^. For simplicity, let 63 = 1, = —i^p and. 



F{d) = dpp' I duju;' 





+ eiP ^2 + e2P ^_2i^pd/c ' ^ 
si - eip' S2 - eip' 



+ , Sl+P S2+P ^_2iujpd/c ' ^ 

Si -p' S2 -p' 
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To evaluate this integral, we could use a; — > as before. But, coth has poles on 
the imaginary axis at Un = "^^^bT ^ ^^^^ integer n. So the closed curve C will 

now have semicircles around these poles in the path that runs along imaginary axis. 
The integration along the semicircles about Un contribute in times the residue of the 
integrand at Un for n > and at n = it will contribute y from the quarter circle. 

h 2TTikBT 



Fid) 



27r2c3 



oo 
n=0 



dp p 



{_einP_ S2n + e2nP ^2ir.pd/c _ ^ 
I Sin - einP' S2n - GlnP' 



^ I Sm + P S2n + P _^2i„pd/c _ ^ 
Sin -P' S2n - p' 



-1 



(A.77) 



where Sjn = a/p^ — 1 + £jn, Cjn = Gjii^n), j = 1,2 and the prime on the summation 
sign indicates that a factor half must be included in n = term. This result was first 
derived by Lifshitz [4J. For the case when ei^2 oo. 



F{d) 



2kBT 



n=0 

oo 



TTC 

2kBT c 



n=0 

.3 °° 



dp p^ 

X) 

dp 



2 



p- 



n=0 



2 



dy 



y 



ey - 1 



where x = ^^^^'^ ■ When x ^ 1, the dominant contribution if from the n 



(A.78) 

(A.79) 
term 



- (A,81) 



Note that the distance dependance changes from \jd^ at 'low temperature' to l/rf^ 
at 'high temperature'. Also, the high and low temperature regime are defined by the 
dimensionless parameter x = '^^^^'^ ■ At any given temperature T, we can always 
choose a large enough ci, such that a; ^ 1. Thus, low or high temperature is equivalent 
to smaller or larger separations. 
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B. CASIMIR FORCE BETWEEN 

SPHERICAL LENS AND FLAT PLATE 



The Casimir force per unit area between two parallel conducting plates, at T = K 
is given by, 



f (z) = (B.l) 

0-013 , _2 , ^ 

= ^ — dyn.cm where Zp, = z va. microns. (ij-2) 

For a finite temperature T we approximate the distance dependence as 

fcT^z) = fcoiz), for z < At; (B.3) 

fcAz) = fcoiz)^, for z>Xt. (B.4) 

AT 

For convenience, we write Eqns. IB.ll and lB.4l as. 

fco{z) = --^ dyn.cm"^ (B.5) 
fcT^z) = dyn.cm"^ (B.6) 

where Ac = — : — . (B.7) 
240 ^ ^ 

The Casimir force between the lens and plate can now be estimated using the "Der- 
jaguin (proximity) approximation", often called "proximity theorem" [T]. Here the 
total force is calculated by integrating over concentric rings on the lens surface that 
are equidistant from the plate. Thus, the Casimir force between the lens and plate is 
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Fig. B.l: Reference axes for the system. 



given by, 



Uz) 

with {2R — x)x 



f{z + x{r)).27rr.dr, 



X 



2R' 



(B.8) 
(B.9) 
(B.IO) 



where z is the separation between the lens and the plate, R is the radius of curvature 
of the lens and rmax is the aperture radius of the lens [Fig. IB.lj . The function, f{z), 
takes different forms depending on z as compared with At- For z < At, 



fcl{z) 



fco{z + x{r)).2nr.dr + 

Jy^mXr 

V^RX¥ ^fi^ rr^a. rdr 

2nAc j+ / — ^ 

[z+Q JV2RX^ At + 

27TRAr ( 1 1 



fcT{z + x{r)).2nr.dr (B.ll) 

fB.12) 



Z-' (Z + At)3 



27rRAr 



2A 



T 



[z + Xry 



z + 



2R 



Uz) 



27lRAr 



iz + Xi 



+ 



2Xi 



[z + XtY {z + X, 



(B.13) 
(B.14) 
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where 



' max 

2R ■ 



For z > At, 



fcliz) 



pr 

Jo 



fcT{z + x{r)).2TTr.dr 



ttRAc 
Xt 



(B.16) 
(B.17) 
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C. YUKAWA FORCE ON A GOLD COATED 
LENS DUE TO GOLD COATED PLATE 



Let the Yukawa potential due to a point mass m at a distance r from the mass be 
given by, 

Y(r) = -niGap . (C.l) 

r 

Consider a glass plate of thickness, g and density, p„, coated with a thin layer of gold 
of thickness, A and density, ps- The Yukawa potential due to the plate at a height, 
z [Fig. IC.lj from the surface of the plate is obtained by integrating the contribution 
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Fig. C.l: Schematic of the plate with coating. 
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due to various mass elements on the plate, 

r g-r/A 

U{z) = -Ga / p(r) —^dV; 
where = w'^ + z'^ . 



U{z) = —Ga / p{z') — =2ttw dw dz 

J a/(w2 + z'^) 

fZ PCX) 

= -Ga2n / dz' p{z') / dr e-"'^ 

Jz+g+A. J z' 

= -2nGaX [ dz' p{z')e-'' 

Jz+q+A 



'z+g+A 

-2TrGa\ (p., [ ^ dz' e~''/^ + Ps [ dz' e""'/^ 

V J Z+g+A Jz+A 

27iGaX^ [p, {e-(^+^)/^ - e-(-+^+^)A} + 



For A << g, 



C.2) 
C.3) 

C.4) 
C.5) 
C.6) 
C.7) 

C.8) 
C.9) 



Uiz) ^ 27rG«A2[p,e-(^+^)/" + p4e-^/"-e-(^+^)/^}]. (C.IO) 

Keeping in mind that our interest is in the region A << we use Eqn. IC.lOl in further 
analysis. The force due to this potential is, 

dU 

m - (c.ii) 

= 27rGaA[p,e-(^+^)/" + p,{e~^/^-e-(^+^)/^}] (C.12) 

= 27cGaX [ps - {ps - Pv)e-^/^] e~'l^ (C.13) 

= /ce-^/^ (C.14) 

where = 2'nGa\ [p, - (p, - pje"^/^] . (C.15) 

The lens has a radius of curvature i?, a volume density p/ and a surface coating of 
uniform thickness b with density pc- We now proceed to derive the force on the lens 
in two limiting cases: (a) A << 5 and (b) A >> b. 



Case (a): \ « b Here the lens may be assumed to be completely made of the coating 
material alone. 
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Fig. C.2: Schematic of the lens with coating. 



Ta{z) = j fiz') pe dVi; (C.16) 

= j fce-^'^^^'^Pc S dAix); (C.17) 

where = 2Rx. (C.18) 

P'^ max 

Ta{z) = fcPcS e-'/' / e-^'/2^"27rrrfr (C.19) 

Jo 

= fcPcS e-^l^ 2'kR\ (l - e-^--/2^^ j (C.20) 

For rlJ2R»X, (C.21) 

^ f,27rRXp,6 e''/\ (C.22) 



Case (b): X >> 6 This case contains two terms one for the contribution from the 
coating and the other for the body of the lens. 

bl. Contribution from the coating: This is same as in case (a) and is given by 
EqnOl 

J-w(z) = 27r RXp, 6 e-'/\ (C.23) 

b2. Contribution from the lens: 

Tm[,z) = J f{z) pi dVi- (C.24) 
= y /, e^(^+^+^')/^ pi nr^dx, (C.25) 
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/, e-^'+^^'^ pi / e-^/^ 2nRxdx (C.26) 



(C.27) 



/e e-(^+^)/^ pi 27rRX 



r.2 



2R 

For rU2i?»A, (C.28) 
J'b2{z) « /e e-(^+^)/^ A 27ri?A. (C.29) 

The net force due to Yukawa potential between the lens and the plate for A >> 5 is 
given by, 

J^b = ^61 + ^62 (C.30) 

= U 27rRX {pc5 + piX e-'l^ ) e-^/\ (C.31) 

In the case our experimental set up, the lens and the plate are both made of the same 
material (BK7 glass) and coated with the gold of equal thickness (1 pm). Hence, 

Pi = Pv ] Pc = Ps and A = 6 (C.32) 

Thus, the force in the two cases is given by, 

J^aiz) = Air'^GaRX^ [ps - (ps - Pv)e~^^^] Ps S e''^^, for A « 5. (C.33) 
J^b{z) = iTT^GaRX^ [ps - {ps - Pv)e-'/^] { Ps 6 + p, X e"^/^ ) e-'/\ {CM) 

for A >> 5. 



